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Abstract 


Benton, James R. Design and Navier-Stokes Analysis of Hypersonic Wind 
Tunnel Nossles. (Under the direction of Dr. John N. Perkins) 

Four hypersonic wind tunnel nossles ranging in Mach number from 6 to 17 are de- 
signed with the method of characteristics and boundary layer approach (MOC/BL) 
and analysed with a Navier-Stokes solver. Limitations of the MOO/BL approach 
when applied to thick high speed boundary layers with non-sero normal pressure 
gradients are investigated. Working gases include ideal air, thermally perfect ni- 
trogen and virial CF^. Agreement between the design conditions and Navier-Stokes 
solutions for ideal air at Mach 6 is good. Thermally perfect nitrogen showed poor 
agreement at Mach 13.5 and Mach 17. Navier-Stokes solutions for CF\ are not 
obtained, but comparison of the effects of low 7 to those of high Mach number 
suggests that the Navier-Stokes solution would not compare well with design. 
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1 Introduction 


The recently renewed interest in hypersonic research, due to such projects as the 
National Aero-Space Plane (NASP) and the Aero-assisted Orbital Transfer Vehicle 
(AOTV), has reemphasized the need for state-of-the-art test facilities in the hyper- 
sonic flight regime. However, a steady decline in funding in this area over the past 
25 years has left the ground testing community unprepared for the challenge. Since 
the early 1970’s, the number of active hypersonic wind tunneb has dropped from 70 
to about 15. Only one major hypersonic wind tunnel has been developed during this 
period and many of the remaining facilities are in need of upgrades [1]. Of prime 
importance in any aerodynamic wind tunnel operation is high quaHiy flow. That is, 
flow in the tunnel test section which is highly uniform with regards to Mach num- 
ber, total pressure and flow angularity in the radial, axial and transverse directions. 
The necessity for uniform flow in ground testing facilities is well understood, but as 
methods of Computational Fluid Dynsunics (CFD) continue to mature, the benefits 
of flow uniformity to CFD grow clearer. Navier-Stokes and Euler codes developed 
for supersonic and hypersonic flow analysis generally hold upstream botmdaries con- 
stant at some uniform freestream condition. So, the ability to deliver uniform flow 
grants the hypersonic wind tunnel the additional objective of validating computer 
solutions. In turn, each time computational data is confirmed by experiment, con- 
fidence in CFD is elevated and the ability to provide advanced designs is enhanced. 
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The fonnidable task of providing high quality flow in a wind tunnel operation is 
achieved almost exclusively through nozzle design and construction. Presently, of 
the seven operational hypersonic wind tunnels at NASA Langley’s Hypersonic Facil- 
ities Complex (HFC), three have unacceptable flow characteristic in the test section 
— the Hypersonic Nitrogen Tunnel, the Hypersonic Helium Tunnel, and the Hy- 
personic CF\ Tunnel. The importance of these tunnels to the maintenance of a 
continuous and well-rounded testing facility at HFC prompted NASA to propose 
upgrades to these and other wind tunnels. Improved flow quality through nozzle 
redesign for three HFC wind tunnels has recdved highest priority in these upgrades. 

In addition to redesigning these nozzles, two new nozzles are proposed for manu- 
facture. The first is a Mach 6 air nozzle for an undeveloped wind tunnel which is to 
be re-machined out of an existing nozzle that was originally designed for Mach 10. 
The second is another nozzle for the Hypersonic Nitrogen Tunnel designed for Mach 
13.5. Design conditions and constraints for all the proposed nozzles are presented 
in Table 1. 

The H3rpersonic Nitrogen Tunnel is an axissrmmetric blowdown tunnel with an 
open-jet test section. The nozzle is 10.5 feet long with a 16 inch diameter exit 
dedgnad for Mach 17 flow. The results of a tunnel calibration, Fig. 1.1, show 
a severe disturbance in the test section [2]. Irregularities are evident up to two 
inches across the centerline, v'hich, in conjunction with the tunnel’s inherently thick 
boundary layer, limits the size of models to about l| to 2 inches in total thickness [3]. 
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At the test section of the CF\ tunnel, a centerline disturbance appears as a Mach 
number spike corresponding roughly to a 13 percent rise in pitot pressure [3]. The 
pitot rake profile of the test section is plotted in Fig. 1.2 and shows uniform flow 
across the asdsymmetric nozzle with the exception of the centerline [5]. Thompson 
and Sutton [6] discovered a discrepancy existing between the reference enthalpy of 
the Method of Characteristics and boundary layer codes used in the nozzle’s design 
by Johnson et al. [7]. This inconsistency produced errors in the density profile of the 
boundary layer which, in turn, gave an erroneous displacement thickness calculation. 
The centerline disturbance, which is the result of the error in displacwnsnt thickness, 
severely limits the physical size of models tested in the facility since any model 
must be situated between the centerline disturbance and the wall boundary layer. 
For the 20 inch diameter nozzle with roughly 1.5 inches of boundary layer at design 
conditions the model size is decreased from 18.5 inches thick to less than 9.25 inches 
thick in the presence of the centerline disturbance. 

All nozzles were designed using the classical approach of iteration between a 
Method of Characteristics (MOC) code and a state-of-the-art boundary layer code 
until aU design constraints were met. This iterative design approach will be re- 
ferred to as MOC/BL. The theory used for the MOC design procedure is based 
on the method published by C. B. Johnson et al. [8] , the details of which are left 
to Appendix A. Figure 1.3 is a schematic of the invisdd portion of a supersonic 
nozzle that shows the various flow regions. The line DE represents the flnal char- 
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acteristic along which the Mach number equals the design Mach number and the 
flow angularity is zero. The region BCD is known as the radial flow region and the 
centerline Mach number here is calculated using radial flow equations in the design 
phase. With the MOC approach, one must specify both the distribution of Mach 
ntimber or velocity on the centerline between points A and B and the maximum 
turning angle at the inflection point C. These parameters effect both the shape and 
length of the expansion region [6]. Boundary Layers were generated with a robust, 
state-of-the-art code developed by E. C. Anderson [9]. Some details concerning the 
boundary layer code used here are presented in Appendix B. 

MOC analysis of the CF 4 and Mach 17 nitrogen characteristic contours were 
performed during the course of this research, but are not included in this thesis. 
For both cases, the MOC analysis procedure agreed with the design. Also, analysis 
of the existing CF4 nozzle at off design conditions was performed and agreement 
with experiment was qualitatively good. Further confidence in the MOC/BL proce- 
dures was obtained by Thompson and Sutton’s [6] prediction of the Mach number 
spike present in the original CF 4 tunnel by an MOC emalysis approach. Also, an 
independent Euler analysis of the Mach 17 Ni characteristic contour gave uniform 
Mach 17.2 flow when it was started with a supersonic inflow profile generated by 
the MOC design code. 

Confirmation of the designs by an MOC zmalysis procedure, however, proves only 
self-consistency between the two MOC/BL design and analysis procedures and not 
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absolute confidence in the designs. And even with certain independent validations 
such as off-design conditions and Euler solutions, there still remains some question 
as to the applicability of the MOC/BL procedures to hypersonic nozzle design with 
its thick, highly turbulent, supersonic boundary layers. The philosophy behind 
MOC/BL as applied to viscous flow centers around the idea of the displacement 
thickness, 6*. Figure 1.4 schematically illustrates the idea behind 6*. In words, 
it is the distance that the actual wall boundary would have to be ^placed such 
that the mass flow in a uniform tnvuctd boundary layer profile, saoving ai the edge 
velocity, is just equal to the mass flow in the oripnal viscous profile. At this new 
displaced wall, the characteristic waves are assumed to reflect or cancel, depending 
upon the local slope of the new invisdd wall. 





— y cos 6) dy — 


Awo 

P.u. 


( 1 . 1 ) 


It is calculated by trapezoidal integration of eq. (1.1) in the boundary layer starting 
at the wall until the integrand is close to zero at the boundary layer edge. This 
technique works well for thin boundary layers because the displacement thickness — 
where invisdd characteristic waves are assumed to be reflected — and where the char- 
acteristic is actually reflected within the viscous boundary layer, are more likely to 
be near the same location. As the boundary layer grows thicker, or as the edge Mach 
number increases, more supersonic flow is found inside the boundary layer and per- 
haps even inside 6*. In the supersonic boundary layer the characteristics would tend 
to curve towards the wall as the Mach number decreased and could therefore reflect 
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from the wall in a different axial and radial location than the 6* approach predicts. 
An error such as this would undermine the relationship between the physics of the 
flow and the mathematical model. For similar reasons the viscous/invisdd approach 
is more accurate for laminar boundary profiles than for turbulent ones since turbu- 
lent boundary layers generally have steeper velocity profiles where the flow could 
remain supersonic very close to the wall. For hypersonic flow, MOC predicts that a 
long, slow tiiming contour is necessary for the proper cancellation of characteristic 
waves to produce uniform flow at the nosale exit. Because of the extended length 
of these nosales, the boimdary layer growth is significant. Another question arises 
here as to the validity of using MOC and boundary layer approach. Specifically, as 
the viscous interaction parameter increases as it does in hypersonic fiow and as the 
boundary layer thickness increases, then the assumption of zero normal pressure 
gradient used in the boundary layer equations becomes less viable. For the present 
designs, the boimdary layers are thick, turbulent and contain a large proportion of 
supersonic flow. They are also characterized by exit Mach numbers between 6 and 
17. Therefore, a Navier-Stokes analysis of the nozzles is desired as an independent 
confirination of the designs. 

Design and Navier-Stokes analysis results of four of the proposed 5 wind tunnel 
nozzles are discussed in this thesis. These nozzles are designed to deliver uniform 
flow of ideal air at Mach 6 (Air) , thermally perfect at Mach 13.5 (ATj — 13) and 
Mach 17 {Ni — 17) and viiial CF^ at Mach 6 (CF4). The objective for this research is 
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twofold. Obviously, bigh quality designs are desired to upgrade HFG test facilities. 
But inherent in this effort is the second objective to survey and evaluate current 
design capabilities as applied to high Mach number contoured nozzles. 

The Navier-Stokes solver, developed by Adjay Kumar [11], was obtained from 
LRC and the appropriate changes were made to accommodate the three gases. Al- 
though the thermodynamic character of the present working gases vary extensively, 
one Navier-Stokes solver was developed that can run these and virtually any single 
specie gas with a bare minimum of changes between cases. This code utilises tabular 
data to replace the respective equation of state throughout. The exact same thermo- 
dynamics that are used in the MOC procedures are used in the Navier-Stokes solver 
in this approach, therefore eliminating errors due to inconsistent thermodsmamics. 
The approach does, however, require the development of a separate computer pro- 
gram or driver for each gas to establish the required tabular equation of state. This 
approach is necessary when analyzing the CF^ nozzle because of its non-standard 
equation of state. CF 4 has a virial equation of state and can be solved explicitly only 
for pressure. To back out other state variables would certainly require an iterative 
procedure or, as in this procedure, table look-up. 
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2 Governing Equations 

The ajdsyininetric Navier-Stokes equations (2.1), are used to describe the flow- 
field. These equations are parabolic in time and elliptic in space and are written 
here in weak conservation form. 

Ut + M, + -{yN)y + - = 0 (2.1) 

V y 

where 

P 

pu 

pv 

pE 

pu 

pu^ + (T„ 
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pUV + T^ 

(pE + <r„)u + T^v + g. 

pv 
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N = 

pV^-\-<Tyy 
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and 


H = 


0 

0 

—att 

0 


E = e + j(tt’ + V*) 


= P+llfe(g-2g + 5) 

^vy - P+|&(£-2g + }) 

<th = P + + 6 “ 2 }) 

= “■ilr(S + fe) 

9- - + 

9w = - (^ + ^) SSS 

fi = fM + fh 

P. = %i 

For the ideal gaa and CF4 case, the above equations are written in terms of the fol- 
lowing non-dimensional variables. For the thermally perfect JV3 cases, diihensional 
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forms of the above equations are used. 


X = xjLc 

y = y/Lc 

u 

II 

p = p/peV^ 

P = pipe 

T = f/Tc 

II 

tS 

II 

h = i/V.’ 

M = A/Me 

Re = 

PcVcLc 


Me 


where the values are dimensional and dimensional constants are defined as, 

Lc = noszle throat radius 
Vc — sti^ation speed-of-sound 
Pc = stagnation density 
Te = stagnation temperature 
Pc = stagnation viscosity 

To complete the set of governing equations, appropriate gas models are needed for 
each type gas — ideal, thermally perfect, real-gas (see Chapter 4) — as well as 
relations for the transport properties, pi and Pr. For ideal gases, pi is calculated 
using Sutherland’s viscosity law and Pr is held constant. For thermall perfect and 
virial gases, the laminar transport properties are interplolated from precalculated 
tables. 

Appendix C illustrates the procedure for determining the turbulent transport 


properties, Pr« and pt. 
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The original Navier-Stokes solver was developed for two dimensions [11] and the 
equations were integrated in strong conservation law form (2.2). 


Ut + M, + Ny = 0 (2.2) 

An axisymmetric patch was later ammended to the code and the governing equa- 
tions are now in weak conservation law form (2.3) Equation (2.3) is obtained by a 
rimple application of the chain rule to the third term of equation (2.1). 

ITt + M. + JV„ -I- * 0 (2.3) 

y 


r 
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3 Integration Technique 

The nozzle flow domain is discretized with a grid similar to Figure 3.0a. An al- 
gebraic numerical transformation provides direct grid control through a stretching 
parameter [11]. Grid points are compressed near the wall to res<flve the boundary 
layer and near the centerline to capture the mathematical discontinuity. Governing 
equations are transformed to a unifomoly spaced grid and integrated using MacCor- 
mmack’s explicit predieior-correetor scheme. This scheme is second order accurate 
in time and space and straightforward to program. In the transformed plane, ((>*7)? 
the governing equations become, 

^ + + = 0 (3.1) 

where the Jacobian of the transformation, J is deflned, 

The unsteady equations (eq. (3.1)) are stepped through time using the two step 
procedure, eqs (3.2) and (3.3), and viscous time step, eqn (3.3), until a steady flow 
condition is achieved and the solution is converged. 

[(Mifc, - Nx^r + (^*< - + JT*] (3.2) 

-dt [{Myr, - Nxr,f^ + {Nxi - My^f^ + } 



The convective terms or outer derivatives (the second and third terms of equation 
(3.1)) for the predictor step are evaluated with forward differences at every odd 
time step and with backward differences at every even time step. These terms 
are evaluated with opposite direction differencing for the respective corrector steps. 
The stress derivatives or inner derivatives (these terms are internal to M and N 
and are defined in equation set (4.1)) are calculated with, differences of opporite 
direction than their respective outer derivatives. Thus, tlm liscona stiesses and 
heat transfer in fiux vectors for a forward predicts st^ will be caknlated with 
backward differencing. This reversal of the differencing is found to improve flow 
symmetry at the reflective centerline [11]. 
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4 Gas Model 

Three gas models are used in the Navier-Stokes solver. The first is the simple 
calorically perfect gas characterized by constant Cp and C, and adherence to the 
standard gas law. This model is used to analyze the Air nozzle and while exper- 
imenting with new calculations, turbulence models, grids, differencing techniques, 
etcetera. The second gas model is that of a thermally perfect gas. Here the gas 
cannot be chemically reacting and intermolecular forces are ignored, so again, it is 
assumed that the gas adheres to the standard gas law. The spedfic heats, enthalpy 
and internal energy depend only upon temperature. 


c = c(T) 


h = h{T) 


C,= 


dr 



The third gas model is that of a virial coefficient gas. Virial gases follow an equation 
of state of the general form. 




( 4 . 1 ) 


where v is the molar specific volume, R is the universal gas constant, and H(T), C(T), 
and ■D(T) are temperature dependant virial coefficients. These coefficients are 
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derived through statistical thennodynaznics to account for intermolecxilar forces. 
Thus, virial gases are, by definition, real gases such that. 


e = e(p,T) 

h = h{p,T) 



4.1 Nitrogen Nozzles 

The iV) — 17 and — 13 are analyzed using the thermally perfect gas model. 
For this case, transport properties, m and Prj, and relationships for h, e and Cp 
are tabulated as functions of temperature. The tables are used at every predictor 
and corrector step during the decomposition of the solution vector, U, using the 
relationship between e, the independent variable and T. The transport properties 
{ire updated every twenty iterations. 

4.2 CFi Nozzle 

The nozzle for the Hypersonic CF 4 Tunnel is analyzed using the real gas rela- 
tionships. Because intermoleculztf forces in CF 4 gas (tetrziflouromethane) remain 
significant throughout its vapor state, the following 13 coefficient virial equation of 
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state is required. 


p __ _l •^2 "f” A3 + B 3 T + Cse 

~V-b (V - 6)2 (V - 6)3 

I B4T ^ At + B$T + ,/ A n m\ aV / 

+ (V - 4)‘ ^ ' (v^ri)! 

The variable b is intended to correct for the volume occupied by the gas molecules. 

CF 4 is used as a wind tunnel gas because it is heavy {M = 88 . 01 ) and allows 
high Reynolds numbers (Re) to be achieved. Real gas relationships for enthalpy 
and entropy are given by equations (4.3) and (4.4). The reference conditions for 
the enthalpy and entropy equations are defined as, 


Pr = 6894.8^ 

Tr = 455.6if 

V, = 1/p, =6.2428^ 

K = 4.6571x10*5^ 

Sr = 3542 . 0 j^ 


I 

I 

I 

I 


( 
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A - fc, » «4 (t - Tr) + y (T* - T*) + i(T* - r») + - r/) + [.-'f**(XT + 1) - + 1)] 

g»(r - Tr) + - e-*y>) Bt{T-Tr) fli(r - TQ + 

(V,- 5 )* (V-J)« (V,-J)» 

+^..v.,T - r.,) - {«r [i. (^) + ^ - ^] 

* e..-") [5^ * 3^ - ^ ^ 


+3(Ai + B$r + c,.-**') [ — L_— + — L_ L_ L_1 

'La(V,-8)> 3(v,-t)* a(V-P a(V-5)*j 


[a(v, -«)*■’■ 4(v, - i)« " a(v-J)» " 4 CvTI) 4 
+6(it, + B,r + c,. *»’) - 4 (v-i)« " 8(v-i)»] } 

^ I V > t 1 

[a(V,-5)» 40^31)4 3(V-i)* 4(V-Q«J 

HA, + B,r)i [••V(«v _ 1 ) _ ,«v.(*v, - 1 )] + r I* la 

(^ - 3^) - - ,3^] 

-*• [io^ - ct] - [(7^ - vT^] ^ T<-'' - ••”•)} <“> 

Subsonic inflow boundary conditions often use an isentiopic condition. The equation 
for entropy for CF 4 gas is equation (4.4) 

S - 5, = 04 ln 1 + 64(7 - T,) + ^(T» - T>) + ^(T» - 3?) 

+K(e-^^ - [-^ + 1 

^ -6 2(V, -6)» 4(V-r-5)*j 

+filn ( - (B, - XC,*-") L_'\ 

\Vr-bJ '\V-b Vr-bJ 

B» - KC,e'*^ 


2 

B, - KC%e~*^ 


L(v-S)» (v,-S)v 

1 1 i 


-B, 


13(V 


4 [(V-S)* (V,-6)^J 

1 ^ 1 1 j 

-S)» 8(V,-S)»J‘" % <» a/ 


(4.4) 


Constants, and b used in equations (4.2), (4.3), and (4.4) are listed in 


r 
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Table 2. Solving these equations for variables other than those for which they are 
written obviously requires an iterative procedure. Analysis experience in the present 
research indicates that the computational expense in iterating these equations at 
every time step is considerable when compared to interpolation &om previously 
generated data. A computer program developed by Hunt and Boney [12] for use 
in data reduction at the Hypersonic CF\ Tuimel is used to generate the tabiilar 


information. 
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5 Boundary Conditions 

Four sets of boundary conditions are needed when solving the Navier-Stokes 
equations in an axisymmetric system — upstream, downstream, wall and centerline. 
An initial condition is also required. Details of these conditions are presented in 
this section. For discussion of the boundary conditions used in the spadal marching 
techniques of the MOC/BL design and analysis, see sections A.1.2, A.1.3 and A.2.1- 
A.2.3. 

5.1 Centerline Boundary Conditions 

l/i^th internal axisymmetric flows, all disturbances are created by the wall and 
propagate along characteristic lines to ultimately impinge upon the centerline. This 
phenomenon emphasizes the importance of a robust centerline boundary condition. 
Of course the centerline boundary condition must be mathematically and physically 
accurate, but with the possibility of the existence of oblique shock waves, it must 
also be niunerically stable. 

The main difficulty in treating the axisymmetric centerline is the singularities 
wristing in the source term of the Navier-Stokes equation set and in the Jacobian of 
the numerical transformation for y = 0. The singularities restrict integration from 
proceeding normally on the centerline. The fourth term of Navier-Stokes equation, 
eq. (2.3), contains the apparent singularity for y = 0. It is obvious, however, that 


F r 
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this singularity cannot physically exist in axisynunetric flows. In fact, it can be 
removed by taking the limit of the entire equation set as y approaches zero which 
includes a simple application of L’Hopital’s ride. However, even with the apparent 
physical singularity removed from the Navier-Stokes equations, integration is still 
restricted on the centerline. Recall that the Jacobian for the transformation is 

^ = 1 / *<»•»)• 

Considering the transformation implied by Figure 3.1, it is obvious that is zero 
for all points; x does not change with respect to the ^ direction. When one remem- 
bers the symmetry condition that is basic to the axisynunetric system, it becomes 
clear that regardless of the distribution of y in the q direction y, is also zero on 
the centerline. Therefore, J is singular and integration cannot take place on the 
centerline in the transformed plane. 

Presently, the centerline boundary condition is treated using the symmetry con- 
dition existing in the original Navier-Stokes solver of reference [11]. The actual 
centerline is not in the domain, but is surrounded on rither side by a grid line 
thus avoiding the singularity of J and the Navier-Stokes equation set. The govern- 
ing equations are integrated on points above the centerline and the properties are 
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reflected to the points below the centerline to exploit symmetry, as follows: 

Pm,l = 

P««,x = 

pEm,\ = pE^^. 

However, solutions using this technique are often characterised by oscillations in 
flow properties at the centerline, espedaUy where high gradients exist such as shock 
waves hitting the centerline or flow through the throat section of a noasle» Although 
the numerical accuracy is not likely to be maintained through a stroi^; shadi hitting 
the centerline, this is acceptable since a nossle with shock waves is of fittls interest 
in this research. 

5.2 Subsonic Inflow Boundary Conditions 

Unlike supersonic inflow boundary conditions, which are generally held constant 
at some fireestream condition, the subsonic inflow boundary condition influences and 
is influenced by the downstream conditions. This feedback relationship is difflcult 
to accurately modd. Characteristic Theory indicates that three properties be spec- 
ified on the subsonic inflow boundary. A fourth property is extrapolated from the 
interior. 

For this study, 5«, and flow angularity are spedfled while P is linearly extrap- 

olated from the interior of the flow. Entropy on the inflow face is held constant at, 
5«, the stagnation entropy. P and 5« are used to solve for the remaining thermo- 
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dynamic variables. The isentropic condition allows the following adiabatic energy 
equation to be used to calculate the V, velocity magnitude. 

The flow angularity and this velocity magnitude specify the u and v velodties. 

5.3 Wall and Supersonic Outflow Boundary Con- 
ditions 

Extrapolation is used for both the wall and outflow boundaries. The supersonic 
outflow boundary depends only upon upstream conditions. All properties across 
the entire outflow boundary are obtained in this way, so there are inaccuracies in 
the boundary layer region near the wall where the flow is no longer supersonic. 

At the wall, the temperature is known and the no slip condition is used, therefore 
extrapolation of p from the interior fixes the thermod3mamic condition along the 
wall. Boundary layer theory, which assumes that the pressure is constant across a 
thin boundary layer, gives support to the use of extrapolation at the wall. 

5.4 Initial Conditions 

A quasi-two*dimensional isentropic expansion with constant mass flow rate is 
used as the starting solution. Flow quantities at the face of the grid are known &om 
the previously analyzed section or are specified in the case of subsonic inflow to be 
the subsonic solution to the standard isentropic flow equations. For the subsonic 
case, the quantities are kept constant across the inflow plane. The massflow rate 
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is calculated for the inflow plane and kept constant at this value for all subsequent 
grid locations. 

With the geometry of the nozzle and hence the area ratio, A/ A* known, Newton’s 
method is applied to equation (5.1) to determine the Mach number at each grid 
location. 



The isentropic flow relationship, eqn. (5.2), is used in conjunction with the reservoir 
temperature to calculate 21*, where t represents axial grid location. 

§ = 1 + (*- 2 ) 

The velocity vector is kept parallel to the axis of symmetry for the initial solution 
with a magnitude calculated from the speed of sound and Mi with eqn. (5.3). 

Ui = Midi = Miy/yRTi (5.3) 

Finally the thermodynamic state is fixed by the massfiow rate equation, (5.4), and 
equation of state, eqn. (5.5). 

m 

~ UiAi 



Pi = piBTi 


(5.5) 
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6 Results and Discussion 

The M0C/6L design approach requires a specified wall temperature and cen- 
terline Mach number distribution in the initial expansion region (between points A 
and B of Fig. 1.3). The wall temperature distributions for the air nozzle and the 
nitrogen nozzles are based on a heat transfer study done on the wall of a Mach 17 
hypersonic wind tunnel nozzle by Sverdrup and are plotted in 6.0. Notice that 
the wall temperatures in Fig. 6.0 start at about lOO^R to 200*12 below the stag- 
nation temperature at the throat and rapidly decrease to some specified constant 
temperature. The decrease follows the behavior of the freestream. The chcnce of 
AB Mach number distribution is generally linear with the exception of the iV] — 17 
nozzle. 

The general approach for the Navier-Stokes analysis of the present designs was 
to divide the nozzle into short sections and analyze them separately. Attempts to 
load an entire nozzle into a single run failed because the CFL stability condition 
near the throat was drastically different than that near the nozzle exit. This was 
due to the la^ gradients in the throat section compared to those near the exit. 
Since the Navier-Stokes code used here has no capacity for adjusting the CFL 
condition at different grid locations the only alternative was to analyze the nozzle 
in shorter sections. Each time a new section of the nozzle was analyzed, the flow 
domain was overlapped to account for the upstream propagation of information in 
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the boundary layer region and to reduce the influence of the extrapolation at the 
outflow boundary on the final solution. In areas of high wall curvature and high axial 
gradients, such as near the throat and inflection point, 0.1 inch x grid resolution was 
used. Downstream of the inflection point where the pressure gradients are low, 0.2 
inch or more x grid spacing was used. The first section (throat section) employed 
the subsonic inflow boundary condition discussed in section 5.2. 

6.1 Mach 6 Air Nozzle 

The Experimental Aerodynamics Branch of NASA Langley is in po s s ess ion of a 
nozsle from a dismantled Mach 10 air wind tunnel. Instead of scraping the nozzle, 
it was decided to use it as the raw casting for the nozzle of an undeveloped Mach 6 
wind tunnel. The new contoured nozzle must fit inside the existing casting which is 
1 inch thick for the majority of its length. The total length is fixed as is the nozzle 
exit diameter, and it must deliver uniform Mach 5.95 flow. 

The contour was designed for perfect gas, 7 = 1.4 since air behaves perfectly at 
and below the stagnation conditions, Po = psia, To = 1260° i2. The wall temper- 
ature distribution was fixed to that shown in Fig. 6.0 and the AB Mach number 
distribution was specified as linear. Several MOC/BL iterations were performed 
where the maximum turning angle, was varied until the contour fit inside the 
tolerance and exhibited the design exit diameter of 14.4 inches. At this point 9mam 
was adjusted to 12 ° so that a TniniTmiTn of metal will be removed from the existing 
nozzle. The total length of the nozzle at this point was 66 inches while the nozzle 
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cactii^ M 76 inches. So, before adding the boundary layer correction for the final 
iteration, a straight, horizontal 10 inch long section was added to the character- 
istic (invisdd) coordinates and the edge conditions were kept constant. The final 
contour is compared to the original contour in Fig. 6.1.1 and is tabulated in Ta- 
ble 3. The section numbers indicate the three nozzle sections and where they bolt 
together. A minimum of 0.0054 inches and a maximum of 0.4524 inches of metal 
will be removed. Fig. 6.1.2 shows the inviscid coordinates and boundary layer edge 
in relation to the final contour. Calculations of 5, based on the velodty profile, 
predicted deven inches of isentropic core at the nozzle exit. 

The Navier-Stokes analysis of the air nozzle was done in 8 sections, each \tith 
0.1 inch longitudinal spacing and 51 points in the radial direction. The final grid 
network was 769 x 51 points for the entire nozzle. Calculated static pressure and 
Mach number contours are shown in Fig. 6.1.3 and Fig. 6.1.4 respectively. The 
Mach 6 contour of Fig. 6.1.4 clearly shows the formation of the desired uiuform 
flow core along the final Mach line corresponding to line DE of Fig. 1.3. The 
horizontal portion of the Mach 6 contour near the nozzle wall shows the straight 
10 inch section added to the inviscid coordinates of the nozzle. As expected, the 
flow in this r^on remained fairly uniform. While the design is for Mach 5.95 
the Navier-Stokes solution predicted Mach 6.06 flow across the nozzle exit plane, 
a 1.85% difference. Corresponding to this. Fig. 6.1.3 shows the static pressure in 
the uniform flow region oscillating about 185 Pa which is 10.5% lower than the 
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design pressure of 207 Pa. A small compression wave is seen in Fig. 6.1.3 starting 
at the wall and impinging upon centerline about halfway down the nozzle. This 
point on the centerline corresponds to point D of Fig. 1.3 and any disturbance 
occurring here can be traced back along the characteristic CD to emanate from the 
inflection point, C. Directly downstream of the compression on the centerline, the 
flow expands upward into the core region causing some oscillations of static pressure 
and Mach number. The downstream effects of this compression/expansion appear 
to damp out as Fig. 6.1.5 shows the exit plane Mach number profile where the 
variation is within 0.25%. Fig. 6.1.6 shows the centerline ratio (P/P«)it 5 compared 
to design. The rise in pressure ratio near x s: 35 inches is evidence the compression 
wave hitting the centerline. Again, the oscillations in {P/P^n8 uc damped to a 
minimum well before the exit. 

Since all total pressure measurements in supersonic wind tunnels are taken be- 
hind behind normal shock waves, it is important to examine the flow quality behind 
a shock. Exit plane data was used as upstream conditions for ideal gas normal shock 
relations to get conditions downstream of a normal shockwave. Fig. 6.1.7 shows the 
profile of the exit plane ratio, (P«)y/P«, where y denotes downstream of a normal 
shock and x is upstream. The variation across the core is within 0.46% and the 
mean value of (P«)y/P» is 4.1% higher than design. The effects of oscillations in 
upstream static pressure is evident in Fig. 6.1.8 showing the downstream total pres- 
sure (P«)y profile across the core with variations up to 1.16%. Again these small 
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variations across the profiles are a result of the compression/expansion activity at 
and beyond the design point, D. The mean value of (Po)y across the nozzle exit 
plane is approximately 6.5% lower than design. This result is consistent since the 
nozzle was predicted to be overexpanded which should cause the normal shock to 
be stronger and yield lower total pressure behind the shockwave. It is the 10.5% 
lower Pa that makes the ratio, (Po)y/Pa, 4.1% higher than design. 

Fig. 6.1.9 shows the MOC/BL design wall pressure compared to the Navier- 
Stokes wall pressure. Recall that the MOC/BL theory assumes zero normal pressure 
gradients in the boundary layer, so the solid line actually represents the behavior 
of static pressure during the isentropic expansion of ideal air from Mach 1 to Mach 
5.95. The pressure is predicted to be about 12.5% lower than design at the noz> 
zle exit. In fact, after about x = 20 inches, the Navier-Stokes wall pressure is 
consistently lower than design. Recall however that the nozzle is predicted to be 
overexpanded by about 10.5% in static pressure. The remaining 2% difference in 
wall pressure at the nozzle exit may be attributed to a non-zero normal pressure 
gradient across the boundary layer, or even a compilation of numerical errors made 
in the design procedure and Navier-Stokes analysis. An indication of the magnitude 
of the numerical error in the Navier-Stokes solution can be seen in Fig. 6.1.10. The 
plot shows a 1.12% increase in massflow rate &om the nozzle inlet to the outflow 
plane. Since this behavior clearly violates the law of conservation of mass, the errors 
must be attributed to numerics and may have influenced the solution. 
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The Navier*Stokes analysis predicted that the Air nozzle contour expands to 
10.5% below design static pressure and 1.85% above Md. In addition, the nozzle 
is characterized by a small oblique compression wave. Oscillations produced by its 
subsequent expansion oil the centerline damp out before the exit. Overall agreement 
between the design and Navier*Stokes analysis of the Mach 6 air nozzle is good. The 
final contour fits inside the prescribed tolerance enabling its manufacture to proceed. 
The design exhibits all of the exit flow requirements within a few percent variance. 

6.2 N 2 — 13 and N 2 — 17 Nozzles 

The Ni — 13 and Nt — 17 nozzles were designed and analyzed tiring the same 
MOG/BL and Navier-Stokes codes used for the Air nozzle. The codes were coupled 
with a set of isentropic expansion zmd thermodynamic property tables for thermally 
perfect nitrogen. Design stagnation conditions for the Nj — 13 and N 2 — 17 nozzles 
are respectively = 4000 ps»o, = 3400* i2 and P« = 4000 psia, To = 3000*i2. 
For these stagnation conditions in the absence of any strong compression or shock 
wave, nitrogen is non- reacting and intermolecular forces are insignificant. 

Both contours were designed for 0mam = 12*. The AB Mach number distribution 
is linear for the Mach 13.5 case and a tlmd order polynomial for the Mach 17 case. 
The third order polynomial was chosen to lengthen the initial expansion region in 
order to maintain vibrational equilibrium. The wall temperature distributions were 
fixed at those illustrated in Fig. 6.0. 

The usual iteration procedure was performed in order to determine the proper 
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displacement thickness to yield the design exit radius. The final nozzle contours are 
illustrated in Figs. 6.2.1 emd 6.2.2 and tabulated in Table 4 and Table 5. Figs. 6.2.1 
and 6.2.2 show the relation of the characteristic coordinates and boundary layer 
edge to the actual nozzle contours. The boundary layer edge profile, calculated 
based on velocity rather than on density, predicted 9.5 inches of isentropic core flow 
at the nozzle exit for the N 2 — 13 nozzle and 8.8 inches for — 17. However, a 
pitot pressure survey across the exit plane flow would not necessarily predict the 
same amount of core flow. Because the pitot tube essentially measures changes in 
entropy across the profile, it will detect the edge of the wall entropy layer whether 
the increase in entropy is due to the velocity gradient or the thermal gradient. 
Thick thermal boundary layers are characteristic of both the nitrogen nozzles. For 
these nozzles, the wall temperature at the exit, 600° i2, is very high compeired to the 
&eestream temperature of equilibrium nitrogen at Mach 13.5 or Mach 17. Under 
the assumption of constant pressure boundary layer profiles, the high temperatures 
at the wall push the mass away from the wall. The larger the temperature gradient, 
the more the density profile is pushed away from the wall. So for the nitrogen 
nozzles, 6 based on this type of density profile will be larger than 6 based on the 
velocity profile such as those of Figs. 6^2.1 and 6.2.2. In fact, consideration of the 
thermal boundary layer in the case of Ni — 13 predicted the isentropic core to be 
almost an inch smaller than when 6 is based on the velocity profile alone. 

Near x — 105 inches on the — 17 contour of Fig. 6.2.2, a small discontinuity is 
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visible. This numerical error was smoothed out using a least squares fit by trial and 
error and is not evident in Fig. 6.2.3 which shows the smoothed iVj — 17 contour 
alongside its first and second derivative. Smoothing of the N 2 — 13 contour was 
unnecessary as the first and second derivatives of the raw design data illustrate in 
Fig. 6.2.4. It is obvious that discontinmties such as that pointed out in Fig. 6.2.2 
could cause oblique shock waves to develop. However, it is not so dear as to how 
sensitive flows such as these are to wall smoothness. The approach here was to 
provide contours that are smooth in the first and second derivative in hopes to 
avoid the development of shockwaves. 

Navier-Stok<» analysis of the JVj — 13 nozde was done in 6 sections. One tenth 
inch axial grid spacing was used for the first 9 inches of the nozzle. The same basic 
grid structure was used in the analysis of the N^ — 17 nozzle. 

Figs. 6.2.5 and 6.2.6 shows the calculated static pressure and Mach number con- 
tours for the N 2 — IZ nozzle. The behavior of pressure contours show evidence of 
a series of oblique compression waves. The small bubble contours on the centerline 
after point D are relative lows in pressure. In between these are repons of compress- 
ing flow due to the compression waves from the upstream wall and impinging upon 
the expanding centerline. The source of the first compression can be traced back to 
just upstream of the inflection point. These compression waves are expanded off the 
centerline to contaminate the downstream flow. In order to darify the behavior of 
the pressure contours near the centerline, the centerline ratio, P/P« is compared to 
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design in Fig. 6.2.7. The alternating compression/expansion is visible. The Navier- 
Stokes solution follows the design closely for only the first 4 to 5 inches. The critical 
inflection point is located at z = 5.1 inches. {P/Po)ns is consistently higher than 
design downstream of the inflection point. 

The Navier-Stokes solution includes implicitly a representation of S*. That is, 
there exists a place within the Navier-Stokes boundary layer repon where the char- 
acteristic waves are reflected back into the mean flow or cancelled. The asaomption 
in the boundary layer corrections was that properties remain constant at the bound- 
ary layer edge values down into the boundary layer until S* is reached. Hence in 
the MOG solution, the slopes of the Mach lines between S and S* are assumed to 
remain constant at each boundary layer profile. In the Navier-Stokes solution, es- 
pecially where the boundary layers get thick, this assumption may not be correct. 
Characteristic waves entering the boundary layer may curve toward the wall as the 
Mach number decreases. Therefore it is possible that the waves in the Navier-Stokes 
solution curved toward the wall between 8 and 6* and reflected or cancelled in a 
location that was further from the wall than assumed in the design. This wotild 
have two consequences. With the two solutions agreeing well near the throat, thus 
fixing the throat area, the larger S*ifs (effective 8* calculated by the Navier-Stokes 
solution) reduces the area ratio for the Navier-Stokes solution thus increasing the 
pressure and reducing the Mach number. This would explain the behavior of the 
centerline pressure ratio in Fig. 6.2.7, which never. reaches the design expansion. 
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The second consequence is the possibility of 6*ss growing large enough to turn the 
mean flow into itself to form the oblique compression waves seen forming near the 
inflection point. 

The above explanation, if accurate, could also explain much of the noted be- 
havior of the Air nozzle. Specifically, the way the nozzle was predicted to expand 
beyond design. The above argument might hold with the exception that the waves 
reflected or cancelled at a point in the boundary layer closer to the wall thus in- 
creasing the effective area ratio. If this scenario is the case, then it indicates a 
serious limitation for which the present MOC/BL design method does not account. 
Unfortunately, the present study is not well suited for the determination of the 
accuracy of this scenario as this study was ori^ally conceived as a development 
project rather than a research project. 

Two contours were designed for Mach 17. The first was designed for a wall 
temperature distribution that qualitatively followed the behavior of the CF4 curve of 
Fig. 6.0 only starting at 3200** i2. Later, wall temperature data from a heat transfer 
study on a Mach 17 nozzle was obtained and used in the design of a second nozzle 
contour. Wall temperature distributions used in the Air and JV3 — 13 designs were 
also based on this data. In both Nz — 17 nozzles the AB Mach number distribution 
followed a third order polynomial. 

The static pressure contours and Mach number contours for the first N^ — 17 
design are in Figs. 6.2.8 and 6.2.9. Again a series of compression waves starting from 
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around the inflection point impinge upon the centerline just beyond the design point 
D. The centerline ratio {P!Po)ns hi Fig. 6.2.10 shows the expansion interrupted 
just beyond the design point by the strong compression. The fact that the nozzle 
expanded to and beyond Mach 17.1 indicates that the source of the compression 
wave is just downstream of the inflection point. The shock wave in this case is 
strong enough to cause the flow to go from the Mach 17.1 to Mach 15 in just 10 
inches. Since the overall objective is to design isentropic noiiles, the Navier-Stokes 
solution was stopped when this shock wave was discovered. 

Dr. A. Kumar of NASA Langley also analyzed this contour and the centerline 
results are ^ven in Fig. 6.2.11. Kumar’s analysis of the same contour used the 
same Navier-Stokes code as the above cases, but with a difierent source for the 
thermodynamics. The results of Kumar’s case and that of flgure 6.2.10 are practi- 
cally identical for the behavior of the centerline pressure ratio. The Mach number 
and pressure contours shown in Figs. 6.2.8 and 6.2.9 are from Kumar’s run since 
the analysis of the first AT] — 17 contour was terminated in the presence of the 
compression wave. 

Analysis results for the second JV) — 17 contour are even further from design than 
for the previous contour. The most significant difference between the two cases is 
that the first contour was predicted to expand to Mach 17.1 and then shock down 
dramatically, but the second case never reached Mach 17. Fig. 6.2.12 shows the 
ratio {P/Po)ns for the second case. Note that the expansion slows down before the 
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design point and almost levels off until an obvious compression wave increases the 
pressure near x = 40 inches. This indicates that the flow was constricted upstream 
of the inflection point producing a series of weak compressions. Now since the mean 
flow is not moving at design conditions the shape of the nozzle wall downstream 
of this constriction — where the nozzle begins to turn the flow — is not correctly 
contoured to caned the characteristic waves. The Mach number is lower than design 
80 the waves have a steeper Mach angle and thus pile up to form a compresdon wave. 

Overall agreement between the design and Navier*Stokez analysis of these two 
nozzles for the Hypersonic Nitrogen tunnel is poor. The analysis of these contours 
reveals a series of compresdon waves emanating near the inflection point, unpingix^ 
upon the centerline, contaminating and increasing the entropy of the mean flow 
region. 

Two characteristics common to all Navier-Stokes analyses in this report are the 
existence of compresdon waves and their apparent source. In the Mach 6 air nozzle 
the wave is very weak so that nonuniformities produced by its expansion from the 
centerline are damped. For the — 13 and — 17 nozzles the wave is stronger 
and the flow is unable to recover from its violent expansion off the centerline. The 
catastrophic downstream effects of the compresdon wave is apparently amplified 
with increased Mach number. The source of the compresdon waves has been traced 
back in all cases to emanate from near the inflection point. This happens because, 
until the inflection point is reached, the nozzle wall is diverging so that it would 
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take large errors in the boundary layer correction to create a strong compression 
wave. But once the inflection point is reached the nozzle wall is already beginning 
to decrease in slope, so it would take a smaller error in 6 * to create a compression 
wave. Therefore, if the cause of these compressions is based in errors or limitations 
in current 6 * calculations for these high Mach number nozzles, any compression 
wave would most likely appear near the inflection point or downstream. The only 
explanation offered thus far as to how the compressions are formed has targeted 
assumptions made in the design. The concept here has been outfined earlier in this 
section and in Chapter 1 and will not be repeated. But there are other possible 
scenarios. The Navier-Stokes analysis, being rather expensive, resolved the bound- 
ary layer with 20 to 25 points at each station compared to the boundary layer code 
used in the design which uses 100 points per station. The low resolution used in the 
Navier-Stokes analysis may not be adequate for internal flows where boundary layer 
thickness is so critical. The resolution used in the x direction was no greater than 
0.1 inch. Fig. 6.2.13 illustrates how low x resolution in the Navier-Stokes analysis 
grid can alter the shape of the nozzle contour and even move the inflection point, 
both of which could aid in the formation of a compression wave. As the original 
scope this study was developmental in nature, there is insufficient information to 
determine which, if any, of these proposals contributed to the compression waves 
and underexpanded conditions of the present nozzles. Further study is needed to 
deterxnine whether these features are introduced by the Navier-Stokes analysis or if 
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they actually represent a limitation of the current design methods. 

6.3 CF 4 Nozzle 

For the design of the CF^ nozzle, the same MOC/BL codes were again used 
only coupled with a more complex set of thermodynamic and ocpansion tables. 
The design stagnation conditions are = 1600 psta, 7« == 1260*P. Since virial 
effects of CF^ are significant throughout the pressure range considered here, the 
standard gas law, which is used implicitly throughout most Nasia-Stokes solvers, 
was replaced by the thermodynamic tables in the computer codes. The tables are 
essentially a set of temperature tables at a range of different pressure settings, 
which allow thermodynamic properties to be calculated as functions of any other 
two thermodynamic properties. The wall temperature distribution was calculated 
using an equation in [ 6 ] which was in turn based on experimental data from the 
existing CF^ nozzle. The centerline Mach number distribution was specified as 
linear. The design proceeded as before with the final contour converging at the 
design exit radius of 10 inches at a length of 98 inches. The contour is shown in 
Fig. 6.3.1. The first and second derivative of the raw design data revealed some 
radical slope changes, but were eliminated by some minor least squares smoothing. 
The first and second derivatives of the smoothed contour are illustrated Fig. 6 . 3 . 2 . 
The final smoothed contour is tabulated in Table 6 . 

The Navier>Stokes analysis of the CF 4 nozzle was unsuccessful. After months of 
modifying the originai code for a virial gas and the table look-up scheme and debug- 
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ging, the resulting Navier-Stokes code was numerically unstable for virtually any 
CFL condition and damping coefRcient. Different centerline and inflow boundaries 
were tried without success. No success was obtained by analyzing the same contour 
under the perfect gas assumption for a gas of 7 = 1.1. The solution diverged in all 
cases tried. 

Design of the Mach 6 Air nozzle of section 6.1 and design of the CF4 contour 
followed identical procedures and used the same computer codes. The two noz- 
zles were designed for the same exit Mach number, wall temperature distribution, 
AB Mach number distribution, and reservoir temperature. The major differences 
between the two cases are and the working gas. 

Thompson and Sutton [6] showed that the MOC/BL approach works well for CF^ 
when they conducted a MOC/BL analysis and preliminary redesign of the existing 
CF^ nozzle in an attempt to characterize the centerline disttirbance of figure 2. They 
conducted the MOC/BL analysis of the nozzle at P, = 2252 psia, P» = 1742 psta and 
P« = 1515 psto. Comparisons of exit plane total pressure ratios computed by the 
MOC/BL procedure to profiles measured experimentally showed good qualitative 
agreement. In areas not in the region of the centerline disturbance, the profiles 
showed good quantitative agreement. Since the MOC/BL philosophy worked well 
for the Thompson study, which modeled CP4 gas in the same reservoir pressure range 
as the presently designed contour, the philosophy should therefore be applicable to 
the present CP4 nozzle design. Furthermore, if the reservoir pressure and working 
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gw of the CF4 nonle design do not undermine the philosophj behind the MOC/BL 
approach and these are the only significant differences between this design and the 
Mach 6 air nozzle design, then the satisfactory results of the air nozzle analysis may 
be useful in inferring confidence in the present CF^ nozzle design. 

Howerer, inherent in the equations of motion are terms inyolving M and 7 where 
low 7 hw similar qualitative effects as high M. For example, in the simple ewe of 
one-dimenrional isentropic fiow of ideal gases the following relations axe valid. 



Applying these equations to Mach 6, 7 = 1.1 (idealization of CFa case) and Mach 
17 , 7 = 1.4 (idealization of — 17 case), 
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The results for the two cases are vtsj close for the total density equation and 
comparable for the total pressure equation. Since these equations are essentially 
a subset of the Navier-Stokes equations, their effects are implicitly included in the 
Navier-Stokes equations. This crudely illustrates the concern for putting overdue 
confidence in the CF^ design by comparisons to the Air nosale case. It is possible 
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that the apparent amplification of the compression wave in the analysis of the pro- 
ceeding nitrogen nozzles could take place in the CF^ nozzle only through a decreased 
7 rather than an increased Mach number. 
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7 Conclusions 

1. The Mach 6 air nozzle is designed to fit all design constraints. The Navier- 
Stokes analysis of the design reveals a weak compression wave and predicts the 
nozzle to expand slightly beyond design. Overall however, the Navier-Stokes 
solution shows good quantitative agreement with the design conditions. 

2. The Mach 13.5 nozzle for the Hypersonic Nitrogen Tunnel is deigned to fit 
all design constraints. The Navier-Stokes analysis shows pom agreement with 
the design flow conditions. The nozzle is predicted to be undereiq^aaded and 
has a series of oblique compression waves contaminating the flow. 

3. Two Mach 17 nozzle contours are designed for the Hypersonic Nitrogen Tunnel 
to fit all design constraints. The Navier-Stokes solutions show poor agreement 
with the design flow conditions. Both contours are predicted to have oblique 
compression waves. One contour expands to design conditions on the cen- 
terline before the compression waves shock the core down and the other is 
predicted to be underexpanded. 

4. The Mach 6 CF 4 nozzle for the Hypersonic CF4 Tunnel is designed to fit 
all design constraints. The Navier-Stokes analysis of the nozzle contour is 
imsuccessful because of numerical instabilities in the computer code. Because 
the eflects of low 7 on nozzle flows are suspected to be sitnilar to those of 
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high Mach number ( the — 13 and — 17 cases), the flow in this nozzle 
may be characterized by a compression wave from the inflection point of such 
strength as to shock down the nozzle. 

5. The source of the compression waves seen in all Navier-Stokes solutions ema> 
nating from near the inflection point of the nozzles is not precisely determined. 
If the apparent ampliflcation of the downstream effects of the waves with Mach 
number can be attributed to a limitation in the current nozzle design meth- 
ods rather than a numerical limitation in the Navier-Stokzz analysiz of hi^ 
speed internal flows, then the present state of MOC/BL nozzle design is not 
adequate to give reliable contoured nozzle designs at higher Mach numbers. 

6. If such a limitations exists then it appears to be in between Mach 8 and Mach 
12. The current methods are shown here to work adequately for Mach 6 and 
in an unpublished study for Mach 8. Problems have been recorded ebewhere 
at Mach 12 and in this study at Mach 13.5 and Mach 17. 

7. Further research is required in the area of high Mach number MOC/BL con- 
toured wind tunnel nozzle design in order to determine whether the apparent 
fimitations of the current techniques are rooted in the design methods or were 
caused by the Navier-Stokes analysis. 
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8 Recommendations 

This project was originally conceived as a development project — specifically, 
to develop hypersonic wind tunnels to meet the standard of high quality flow. Over 
the course of the effort many areas within the nozzle design discipline have been 
discovered where limitations in the current methods are apparent. Some serious 
and basic research is needed in this area. The following is a partial list as to where 
the author feels this research should be concentrated. 

1 . Develop a Navier-Stokes solver specifically designed for internal flow problems, 
particularly axisymmetric geometries. This code should have the following 
characteristics: 

• adjustable CFL condition according to gradients within the flow so that 
entire nozzles can be loaded into a single run. 

• use of effective 7*8 and compressibility factor, Z, so that equilibrium 
effects and real gas virial effects can be easily included and updated 
every 5 to 10 time steps. 

• use of finite volume scheme instead of finite difference scheme. The 
author feels that the volumetric approach is a more sophisticated way to 
approach the axisymmetric centerline problem. 
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• research is needed to provide a reliable, accurate subsonic inflow bound- 
ary condition or starting solution for these and other flows. 

2. Perform extensive parametric studies using MOC/BL nozzle design and Navier- 
Stokes analysis to gain basic understanding of how design parameters ef- 
fect nozzle flow and to get a better understanding of the limitations of the 
MOC/BL approach. 

3. If further research concludes that the current MOC/BL techniques become 
inaccurate at some limiting Mach number, research should proceed in devel- 
oping an alternative design algorithm. A spadal marching technique that can 
be used in conjunction with the Parabolized Navier-Stokes equations should 
be investigated. 
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A Method of Characteristics 
Applied to Nozzles 


The equations governing axisymmetiic, irrotational, steady, supersonic flow- 
fields, equations (A.l) and (A.2), are hyperbolic, quasi-linear partial differential 
equations and can therefore be solved using the Method of Characteristics (MOC). 
In a hyperbolic system, there exist characteristic curves (characteristics) along 
which the governing partial differential equations can be manipulated into ordi- 
nary differential equations called eompaU^Uy equations. Across characteristics, 
flow properties may have discontinuous derivatives while the property itsdf remains 
continuous. For these equations. 




tty - = 0 



(A.1) 

(A.2) 


which are a composite of the continuity equation and Euler’s equation (A.l) and 
the irTotaiional condition eqn (A.2), the compatability equations are shown in [10] 
to be, 


(i) 


(A.3) 


The equations are valid along the characteristic lines. 


(^) =t»n(<±M) 


(A.4) 
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Equations (A.3) are readily expressed in finite difference form and can therefore be 
easily integrated over a mesh of characteristic curves described by equations (A.3) 
and (A.4). 

Four equations are available relating the five variables x,y, V, 9 and /x, so one 
additional relationship is needed. This is found by using the definition of the Mach 
angle, 

(if) - (A-S) 

and a real-gas relationship between M and V. The latter is achieved by tabulating 
M and V through an isentropie expansion including real-gas effects. In general, the 
four equations are integrated in a marching scheme using the modified Euler pre- 
dictor/corrector method with an iterative algorithm on the corrector step. Details 
of this procedure are found in [10]. A direct marching technique is used for the de- 
sign and analysis procedures, meaning that the left and right running characteristic 
waves {C+ and (7- characteristics, respectively) make-up the computational grid. 

A.l Design Procedure 

The MOC design procedure used here is based primarily on that of reference 
[8], which is a modification of the work presented in [7|. Integration begins with 
the design condition. That is, the starting solution is the final characteristic DE 
in Fig. 1.3. For this approach one must specify the solution along ABDE. The 
maximum turning angle of a contoured nozzle occurs at point C of Fig. 1.3, the 
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inflection point. The angle of the characteristic wall here, must be specified. 

A. 1.1 Radial Flow Region 


Region II of Fig. 1.3 is called the radial flow region and is governed by the 
characteristic radial flow equations, (A.6) and (A.7), which are the Prandtl-Meyer 
expansion function and the radial velocity relationship. 


and 


Jo ' 2 JVumX V 


, /M 1 yv ^ 

\TcrJ 2Jvu.i V 


(A.6) 


(A.V) 


Here 9f is the flow angle and $t is the Prandtl-Meyer expansion function for radial 
flow. For these equations, r is the distance from the source shown on the centerline 
upstream of point A in Fig. 1.3, and r^. is the sonic radius in source flow. These 
equations are numerically integrated using an arbitrarily small velocity step AV 
and the tabulated relation for M and V until the known design velocity is reached. 
The integration of eq. (A.7) to Vd fixes the location of point D. The relationship 
between V, $i is retained and used as the centerline solution on BD. Using the 
relations, 


9f = — Bt along CD 


(A.8) 


and 


Bf ^ Bt — Big along BC 


(A.9) 
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equation (A.IO) is derived. 

^Ib = ^Id “ 2dc. (A.IO) 

9jg can now be used in conjunction with tbe previously retained relationship be- 
tween $t and ^ to determine the location and flow solution at point B. The cen- 
terline solution is now fixed on AB by assuming a M or F distribution. Using sonic 
conditions at A and the known solution at B, one has the freedom to spedfy a first, 
second or third order distribution. With the solution now fixed along ABD, it is 
left to determine the solution at the exit, the starting point of integration, line DE. 

A.1.2 Uniform Flow Region 

Rerion IV is the uniform flow rerion of a contoured nosrie since it is desired 

I 

to have constant properties and no flow angularity. Therefore, the solution along ' 

DE is known to be the design conditions. However, it is still left to specially locate 
point E, which is fixed by the mass flux through the nozzle since all flow properties 
on DE are known. The calculation of the limiting mass flux employs the stream 
function concept and is reviewed in the following section. 

A.1.3 Stream Function and Inviscid Wall Boundary 

Integration along the characteristics for the derign procedure begins at the cen- 
terline where there are two points of known solution. A C+ characteristic through 
a centerline point and a C- characteristic through a point on the starting charac- 

j 

teristic are used to simultaneously integrate equations (A.3) and (A.4) for the flow 
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solution at another point on the second right running characteristic. The general 
procedure is to calculate the stream function, along with the other flow V2triables, 
at each point until it reaches the vedue limited by the mass flow, is calculated 
in the radial flow region at point C, the inflection point of the contour. Integration 
of the axisymmetric radial flow relation for eq.(A.ll), 

W (A.11) 


where 


? = 


ptUlT%, 


gives the following result after evaluating a constant of integration. 




(A.12) 


Again, equations (A.8) and (A.9) are used with the retained radial flows relationship 
to determine ~ and V at point C. The real-gas relationship between p and V is 
obtained from tables. The equation for ^ used while integrating the characteristic 
mesh is 


Jo Pt 


Fsin^ 



ds. 


(A.13) 


where s refers to the distance along the Mach wave and y is the perpendicular 
distance from the x axis. When, during integration of eqs. (A.13), is reached, 
integration stops, the wall point location is recorded and a new Mach line is started. 
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A. 1.4 Iterative Procedure 

For the present designs the exit radius is a design constraint. This means that 
two nested iterative procedures will be needed to get a converged design. The outer 
iteration is set up to converge on the design exit radius. The inner iteration is 
designed to converge on the correct value of 6* at each x station for each contour 
in the outer iteration. 

An initial guess for 6* at the exit is necessary and the nondimennoxud i&visdd 
contour is scaled so that its exit radius is equal to the design radios minus the 
guessed S*. The boundary layer code (see Appendix B) is used to generate S* 

'C 

for these coordinates and this is added to the scaled invisdd coordinates to get 
an initial guess for the corrected wall, yeorr. Now, each time the coordinates are 
scaled or adjusted the curvilinear distance, s, of the nozzle contour is changed and 
since S* = ^*(s), the displacement thickness is also changed. So, a second iterative 
procedure is followed to find the correct value of S* at each x-location. Here yeorr is 
adjusted at each x-location by 

^wa ll — 2 “ yeorr "I" ^ ) 

after every iteration on 8* until AwaU is less than a specified small e. 

If the exit radius of the newly adjusted wall is not equal to the design radius, 
then this wall is scaled to the design exit radius and the inner iteration loop is 
repeated until they are equal. 
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A. 2 Ancdysis Procedure 

The present MOC analysis code is based on chapters 12 through 15 of [10]. 
Integration of eqs (A.3) and (A.4) along the characteristics proceeds as in the design 
case. The main differences between the design and analysis procedures are in the 
starting solution and the techniques for dealing with the boundaries. Integration 
begins near the throat on what is referred to as the initial value solution. 

A.2.1 WaU Point 

Sketch A.l Hlustrates the wall point situation. The presence of the solid bound- 
ary clearly eliminates the C_ characteristic. However, the wall location is known 
and the flow tangency relationships, B are known (the slope of the wall). So for the 
three unknowns, Af , and V we have three equations, the compatability equation 

for the (7+ characteristic, the tabular real-gas isentropic expansion of the gas and 


the Mach wave relation. 
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= d6 tan p + 


sin p tan p sin 9 
sin(d + p) 



p = 



^ = tanO = ^, 
dx u 


M = M{V). 


A.2.2 Centerline Point 


A witwilaf dtuaiion arises at the axis (see sketch A.2). However for this case both 
characteristics exist. The flow properties at point 1 are merely reflected to point 
2 exploiting the symmetry condition of the axisymmetric centerline. Fiirthermore, 
compatability equations are simplified greatly since y at v are zero at the centerline 
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so that only the characteristic is needed. 



Sketch A. 2 Centerline point solution by MOC. 

A.2.3 Iterative Procedure 

To start the Eoialysis iteration, an initial guess for S* is needed. This is obtained 
by generating a boundary layer using one-dimensional isentropic expansion to the 
design Mach number as the edge conditions. The boundary layer correction is 
subtracted from the actual nozzle wall to get a first guess for the characteristic 
coordinates. The MOC analysis code is run on these coordinates. From the analysis 
output, a second, hopefully better, estimate of 6* is generated. This is iterated until 
the changes in S* are small. 

When and if a strong oblique shock wave forms in the MOC analysis solution, 
then convergence may not be reached. The current status of the MOC/BL anal- 
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ysis procedure is not generally sophisticated enough to account for strong shock- 
wave/boundary layer interactions. 
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2 Numerical Solution of the 
Boundary Layer Equations 

E. C. Anderson developed a computer program with the objective of generating 
turbulent bovmdary layer solutions for hypersonic flows, external and internal [9]. 
A brief overview of the techniques used is presented here. 

The Crank-Nicolson semi-implidt finite difference method is employed to inte- 
grate the steady momentum, energy, species, and continuity equations* The pro- 
gram is limited to gases in chemical equilibrium. By assuming a constant gas com- 
position across the boundary layer, the effects of chemical reactions are included in 
the heat transfer terms of the energy equation. In terms of non-dimensional vari- 
ables these equations are, in axisymmetric form. 


Continuity 

{dpur) d[(pp+^r] 
dx dy 

Momentum 


( 2 . 1 ) 


dn , — du du, d ( du , 

(»^ + (/»>+ 7 ?)^ + ( 2 . 2 ) 


Energy 


dH , — ^dH d 


** -l>VS' + II 


Pr dy 


dy. 


( 2 . 3 ) 
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The aboye non-dimensional, turbulent boundary layer equations are transformed 
using the Levy-Lees transformation (eqs. B.4 and B.5) to appear in the form of 
equations (B.6), (B.7) and (B.8). 



C(») = / p*u^p,r^dx 

JO 

(2.4) 


, ^ p,u,r du, pf P J 

= (2{)* /. 

(2.5) 

Continuity 

2(F( -h F' + 2?* = 0 

(2.6) 

Momentum 

2^FFi + (C(l + e+)f^y 

(2.7) 

Energy 






+^{C7i [C’FF’ + C{F^ + F')] + CFrC[} 

(2.8) 


where, 


V = 


2^ 




u 


F = f = 


. (pv+'p^r 

g=—-, (7 = -^ 
H, * p.u. 


/? = 


2^ du« 


Since Cranlc-Nicolson technique is semi-implidt, the finite difference represen- 
tation of these equations contains independent variables at the station where the 
solution is desired. This requires a matrix inversion, or an iterative procedure. An 
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iterative procedure is used here. The coefficients of the independent variables eire 
evaluated using properties at the previous station. The finite difference equations 
for energy and momentum are solved for the independent variables at the new sta- 
tion. At each iteration the continuity equation is integrated using the trapezoidal 
rule for the normal velocity, v. From the new properties, better estimates for the 
coefficients are calculated and the procedure is iterated until converged. 

Once the solution has converged at a particular station, eq (B.9) is integrated 
trapezoidally from the wall to the point where pujpjj^ = 0.995. 


The point in the boimdary layer where this occurs is called the displacement thick- 
ness and is used as the boundary layer correction to the MOC invisdd flow region. 

The turbulence model is algebraic and is based on a two layer eddy viscosity 
assumption. Hence, the Reynolds stress and turbulent heat transfer are defined as 
follows. 


-7-7 




( 2 . 10 ) 

( 2 . 11 ) 


where. 


pTt = 


Sk 


The total shear stress and heat transfer become. 




( 2 . 12 ) 
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and, 



(2.13) 


The eddy viscosity is modeled assuming two layers. Both iimer layer and outer 


layer models for e'*' are evaluated until 


£.• 



Below this point in the boundary layer ef is used and above this p<»nt e* is used. 

The outer layer is modeled by a zero equation model developed by Clauser (eq. 
B.14). 

/Ml- At 

(2.14) 


- t. <”• 


Here, Sk is the incompressible displacement thickness, €vd u a coordinate stretching 
parameter, and is an empirical constant. 

The inner layer is based on Prandtl’s mixing length theory and uses the Cebeci- 
Smith zero equation turbulence model. The eddy viscosity for this layer is. 


ef = *1 < 


1 — exp 


— *a ^ 


T ■ 


) 


in 


(2.15) 


where, 


jp jkiypu^yr 

cvdM(2^)^ 


= 


yp 


26 fie\n) 


k\ = constant. 
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3 Turbulence Modeling for 
Navier-Stokes Solutions 


The turbulent Navier>Stokes equations are identical to the laminar equations 
with the exception of additional terms in the turbulent expressions for heat trzmsfer 
(3.1) and the stress tensor (3.2). 




'dXm 


I t 

-pUfUj 


(3.1) 

(3.2) 


However, these two terms introduce many additional unknowns to the problem 
which can be defined only by the already unknown turbulent structure. These un- 
knowns are dependant not only upon the physical properties of the particular gas, 
but also upon local conditions such as velocity, geometry, surface roughness, and 
previous spadal and temporal conditions [ 13 ]. Following the original idea of Boussi- 
nesq ( 1877 ), one can model the effect of turbulence — rather than model the actual 
physics — and approach mathematical closure by considering the turbulent terms 
as functions of local geometry and fiow properties. By defining the turbulent shear 
(3.3) and turbulent heat transfer (3.4) in terms of effective turbulent coefficients of 
viscosity, /it and thermal conductivity, ibt. 


K dUi dUj\ 2 dUyn 


(3.3) 


- hc,u\r = 


(3.4) 


64 


the assumption is that the turbulent terms behave similarly to th^ laminar (i.e. , 
physical) counterparts. Thus, in the stress and heat transfer equations of equation 
set (2.1), 

11 = Hi+ fit 


and 


k ^ ki + ki 


In this research the turbulent conductivity is suppfied throi^ the turbulent Prandtl 
number, 

Pt - 

and is considered to be constant at an empirical value of 0.9. The turbulent viscosity 
is furnished by the Baldwin-Lomax two layer eddy viscosity model [14]. 


3.1 Baldwin-Lomax Turbulence Model 


The Baldwin-Lomax two layer eddy viscosity modri was developed to address 
the problem of including turbulent effects in Navier-Stokes simulations of separated 
flows where the boundary layer thickness is not generally known. Baldwin and 
Lomax [14] defined the eddy viscosity on an inner and outer layer by 


P* = 




y ^ Vcrottovar 


(3.5) 


^ (/^)0ulcr y ^ yero§$ov€r 

where yg-oijcw is the shortest normal distance from the wall at which the inner and 
outer coefficients are equal. 


r 


65 


The Prandtl-Van Driest formulation is used for the inner region. 


i"i 


(3.6) 


where 


/ = [l - exp (-i/+/A+)] 


|(i;| is the magnitude of the yortidty 


and 


. . du dv 




The formulation for the outer region is 


(3.7) 


(3.8) 


(3.9) 


{lh)ouUr — KCcPpF^a!mFKlA{y) (3.10) 

where K is the Clauser constant Ccp is &u additional constant and 

yrf%amFmam 

^v$akm ~ ^ ^ (3.11) 

^wkymam'^\if/ ^mam 

^ mm 

The TnaarimiiTn point of the function (3.12) and the point at which it occurs define 
Fn^ and y^. 

F{y) = y Iwj [l - exp (-y+/i4+)] (3.12) 

Fxuk is the Klebanoff intermit tency factor expressed as 


FKUh{y) = 


l+5.5f^)' 

\ ymam / 


-1 


(3.13) 


r 
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and the quantity u^f is the difference between the maximum and minimum total 
velocities, which for no slip boundary conditions is the maximum total velocity in 
the profile. The following constants are used in the turbulent viscosity calculations. 

>1+ = 26 
Ccp — 1*6 
Cxuk = 0.3 
Cwk — 0.25 
k = 0A 
K = 0.0168 
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Tables 


Table 1. Design Conditions and Constraints 


Nozzle 

Working 

Gas 

Exit 

Diam. 

Exit 

Mach 

Total 

Length 

Po 

psia 

To 

N^-n 

Nitrogen 


17 

— 

4000 

3400 

Ni-lZ 

Nitrogen 

20 in. 

13.5 

— 

tfliiinM 

3000 

CFa 

CF 4 

20 in. 

5.95 

— 

1600 

1260 


Heliiun 

21 in. 

20 

124 in. 

1500 


Air 

Air 

~ 7.2 in. 

5.95 

76 in. 

45 

1260 1 


Table 2» Constants for CF^ Thennod 3 m*<T'wc 


R 

94.47002 

FI 

0.0015 

Jk. 

-58.11999 


7.387702 X 10“^ 


0.1032692 

Bz 

3.873415 X 10-^ 

Ct 

-508.9595 

c. 

0.9054616 


2.011774 X 10-« 

At 

-2.929506 X lO"* 

^4 

-7.385870 X 10-* 

Bt 

1.066416x 10-^® 


4.025724x10“ 

Ct 

-3.162000 X 10-^ 

Be 

-1.149705 X 10» 

K 

1.758236 

II ^ 

-1.059139 X W 

- 


C onstants for Enthalpy anc 

1 Entropy Equations 

J 

br temperature interval 100 K to 615 K, 

04 

1.43308 X 10> 

b4 

1.51168 

1 <4 

5.15767 X 10-“ 

d4 

-1.59381 X 10-® 


for temperature inttwal 615 K to 1500 K. 


04 

2.95790 X W 1 64 

1.53754 1 

C 4 

-1.0727 X 10-* 1 <U 

2.64172 X 10-^ 1 
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Table 3. Design Wall Coordinates for the Air Nozzle 


X , in . 

r , in . 

0.0000 

0.9015 

0.1000 

0.9033 

0.2000 

0.9051 

0.3000 

0.9077 

0.4000 

0.9111 

0.5000 

0.9152 

0.6000 

0.9200 

0.7000 

0.9254 

0.8000 

0.9315 

0.9000 

0.9382 

1.0000 

0.9454 

1.1000 

0.9533 

1.2000 

0.9618 

1.3000 

0.9708 

1.4000 

0.9804 

1.5000 

0.9906 

1.6000 

1.0013 

1.7000 

1.0126 

1.8000 

1.0245 

1.9000 

1.0369 

2.0000 

1.0499 

2.1000 

1.0635 

2.2000 

1.0776 

2.3000 

1.0923 

2.4000 

1.1075 

2.5000 

1.1231 

2.6000 

1.1392 

2.7000 

1.1558 

2.8000 

1.1727 

2.9000 

1.1900 

3.0000 

1.2076 

3.1000 

1.2255 

3.2000 

1.2438 

3.3000 

1.2623 

3.4000 

1.2810 

3.5000 

1.3001 

3.6000 

1.3193 

3.7000 

1.3387 

3.8000 

1.3583 

3.9000 

1.3781 

4.0000 

1.3981 

4.1000 

1.4182 

4.2000 

1.4385 

4.3000 

1.4589 

4.4000 

1.4795 

4.5000 

1.5001 

4.6000 

1.5209 

4.7000 

1.5418 

4.8000 

1.5627 

4.9000 

1.5838 

5.0000 

1.6049 

5.1000 

1.6262 

5.2000 

1.6474 


X , in . 

r , in . 

" 5.3000 

1.6688 

5.4000 

1.6903 

5.5000 

1.7117 

5.6000 

1.7333 

5.7000 

1.7549 

5.8000 

1.7765 

5.9000 

1.7982 

6.0000 

1.8200 

6.1000 

1.8417 

6.2000 

1.8636 

6.3000 

1.8854 

6.4000 

1.9073 

6.5000 

1 . 9292 . 

6.6000 

1 . 9 S 11 

6.7000 

1.9731 

6.8000 

1.9951 

6.9000 

2.0171 

7.0000 

2.0392 

7.1000 

2.0613 

7.2000 

2.0833 

7.3000 

2.1055 

7.4000 

2.1276 

7.5000 

2.1497 

7.6000 

2.1719 

7.7000 

2.1941 

7.8000 

2.2162 

7.9000 

2.2385 

8.0000 

2.2607 

8.2000 

2.3051 

8.4000 

2.3496 

8.6000 

2.3942 

8.8000 

2.4388 

9.0000 

2.4834 

9.2000 

2.5281 . 

9.4000 

2.5727 

9.6000 

2.6170 

9.8000 

2.6611 

10.0000 

2.7049 

10.2000 

2.7484 

10.4000 

2.7916 

10.6000 

2.8345 

10.8000 

2.8771 

11.0000 

2.9193 

11.2000 

2.9611 

11.4000 

3.0027 

11.6000 

3.0439 

11.8000 

3.0848 

12.0000 

3.1253 

12.2000 

3.1655 

12.4000 

3.2054 

12.6000 

3.2449 

12.8000 

3.2840 

13.0000 

3.3229 





Table 3. (continued) 


X , in . 

r , in . 

13.2000 

3.3614 

13.4000 

3.3996 

13.6000 

3.4375 

13.8000 

3.4750 

14.0000 

3.5122 

14.2000 

3.5491 

14.4000 

3.5857 

14.6000 

3.6219 

14.8000 

3.6579 

15.0000 

3.6935 

15.2000 

3.7288 

15.4000 

3.7638 

15.6000 

3.7985 

15.8000 

3.8330 

16.0000 

3.8671 

16.2000 

3.9009 

16.4000 

3.9345 

16.6000 

3.9677 

16.8000 

4.0007 

17.0000 

4.0333 

17.2000 

4.0657 

17.4000 

4.0978 

17.6000 

4.1297 

17.8000 

4.1612 

18.0000 

4.1925 

18.2000 

4.2235 

18.4000 

4.2543 

18.6000 

4.2847 

18.8000 

4.3150 

19.0000 

4.3449 

19.2000 

4.3746 

19.4000 

4.4040 

19.6000 

4.4332 

19.8000 

4.4622 

20.0000 

4.4909 

20.2000 

4.5193 

20.4000 

4 . 547 $ 

20.6000 

4.5755 

20.8000 

4.6032 

21.0000 

4.6307 

21.2000 

4.6580 

21.4000 

4.6850 

21.6000 

4.7118 

21.8000 

4.7383 

22.0000 

4.7646 

22.2000 

4.7907 

22.4000 

4.8166 

22.6000 

4.8422 

22.8000 

4.8677 

23.0000 

4.8929 

23.5000 

4.9549 

24.0000 

5.0157 

24.5000 

5.0752 

25.0000 

5.1334 

25.5000 

5.1904 

26.0000 

5.2461 

26.5000 

5.3007 

27.0000 

5.3541 

27.5000 

5.4063 


X , in . 

r , in . 

28.0000 

5.4574 

28.5000 

5.5075 

29.0000 

5.5564 

29.5000 

5.6043 

30.0000 

5.6511 

30.5000 

5.6970 

31.0000 

5.7418 

31.5000 

5.7857 

32.0000 

5.8285 

32.5000 

5.8705 

33.0000 

5.9115 

33.5000 

5.9516 

34.0000 

5.9908 

34.5000 

6.0291 

35.0000 

6.0680 

35.5000 

6.1033 

36.0000 

6 . 1391 ^ 

36.5000 

6.1740 

37.0000 

6.2082 

37.5000 

6.2418 

38.0000 

6.2742 

38.5000 

6.3080 

39.0000 

6.3371 

39.5000 

6.3674 

40.0000 

6.3971 

40.5000 

6.4260 

41.0000 

6.4542 

41.5000 

6.4817 

42.0000 

6.5085 

42.5000 

6.5347 

43.0000 

6.5602 

43.5000 

6.5851 

44.0000 

6.6094 

44.5000 

6.6330 

45.0000 

6.6560 

45.5000 

6.6784 

46.0000 

6.7002 

46.5000 

6.7215 

47.0000 

6.7421 

47.5000 

6.7622 

48.0000 

6.7818 

48.5000 

6.8008 

49.0000 

6.8193 

49.5000 

6.8372 

50.0000 

6.8547 

51.0000 

6.8880 

52.0000 

6.9194 

53.0000 

6.9489 

54.0000 

6.9766 

55.0000 

7.0026 

56.0000 

7.0268 

57.0000 

7.0494 

58.0000 

7.0705 

59.0000 

7.0901 

60.0000 

7.1082 

61.0000 

7.1249 

62.0000 

7.1403 

63.0000 

7.1544 

64.0000 

7.1672 
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Table 4. Design Wall Coordinates for the N] — 13.5 Nozzle 


1 

r, in. 

I 0.0000 

0.1194 

0.1000 

0.1202 

0.2000 

0.1224 

1 0.3000 

0.1256 

1 0.4000 

0.1298 

I 0.5000 

0.1351 

1 0.6000 

0.1413 

1 0.7000 

0.1484 

0.8000 

0.1566 

0.9000 

0 . 1657 . 

1 1.0000 

0.1757 

1.1000 

0.1869 

1.2000 

0.1990 

1.3000 

0.2122 

1 1.4000 

0.2266 

1.5000 

0.2421 

1.6000 

0.2587 

I 1.7000 

0.2762 

1 1.8000 

0.2946 

1 1.9000 

0.3135 

1 2.0000 

0.3331 

2.1000 

0.3531 

2.2000 

0.3735 

1 2.3000 

0.3943 

2.4000 

0.4153 

1 2.5000 

0.4365 

1 2.6000 

0.4579 

2.7000 

0.4796 

2.8000 

0.5013 

2.9000 

0.5232 

3.0000 

0.5452 

3.1000 

0.5673 

3.2000 

0.5895 

3.3000 

0.6117 

3.4000 

0.6341 

1 3.5000 

0.6564 

1 3.6000 

0.6789 

1 3.7000 

0.7014 

1 3.8000 

0.7239 

1 3.9000 

0.7465 

1 4.0000 

0.7692 

1 4.1000 

0.7918 

1 4.2000 

0.8146 

4.3000 

0.8373 

4.4000 

0.8601 

4.5000 

0.8829 

4.6000 

0.9058 

4.7000 

0.9287 

1 4.8000 

0.9516 

4.9000 

0.9745 

5.0000 

0.9974 

1 5.1000 

1.0204 

5.2000 

1 . 0433 L 


X, in. 

r, in. 

5.3000 

1.0661 

5.4000 

1.0888 

5.5000 

1.1114 

5.6000 

1.1339 

5.7000 

1.1503 

5.8000 

1 . 178 S . 

5.9000 

1.2007 

6.0000 

1.2220 

6.1000 

1.2447 

6.2000 

1.2600 

6.3000 

1*2002 

6.4000 

1 . 301 V 

6.5000 

1.3315 

6.6000 

1 . 3 S 2 r . 

6.7000 

1.3742 

6.8000 

1.3954 

6.9000 

1.4165 

7. 0000 

1.4375 

7.1000 

1.4584 

7.2000 

1.4792 

7.3000 

1.4999 

7.4000 

1.5205 

7.5000 

1.5410 

7.6000 

1.5614 

7.7000 

1.5817 

7.8000 

1.6019 

7.9000 

1.6220 

8.0000 

1.6420 

8.1000 

1.6619 

8.2000 

1.6818 

8.3000 

1.7015 

8.4000 

1.7211 

8.5000 

1.7406 

8.6000 

1.7601 

8.7000 

1.7794 

8.8000 

1.7987 

8.9000 

1.8179 

9.0000 

1.8370 

9.1000 

1.8560 

9.2000 

1.8749 

9.3000 

1.8937 

9.4000 

1.9125 

9.5000 

1.9311 

9.6000 

1.9497 

9.7000 

1.9682 

9.8000 

1.9867 

9.9000 

2.0050 

10.0000 

2.0233 

10.1000 

2.0414 

10.2000 

2.0595 

10.3000 

2.0776 

10.4000 

2.0955 

10.5000 

2.1134 



Table 4. (continued) 


X , in . 

r , in . 


X , in . 

r , in . 

10.6000 

2.1312 


21.5000 

3.7423 

10.7000 

2.1489 


21.8000 

3.7796 

10.8000 

2.1666 


22.1000 

3.8166 

10.9000 

2.1841 


22.4000 

3.8533 

11.0000 

2.2017 


22.7000 

3.8898 

11.1000 

2.2191 


23.0000 

3.9259 

11.2000 

2.2365 


23.3000 

3.9618 

11.3000 

2.2538 


23.6000 

3.9974 

11.4000 

2.2710 


23.9000 

4.0327 

11.5000 

2.2882 


24.2000 

4.0678 

11.6000 

2.3052 


24.5000 

4.1026 

11.7000 

2.3223 


24.8000 

4.1371 

11.8000 

2.3392 


25.1000 

4.1714 

11.9000 

2.3561 


25.4000 

4.2054 

12.0000 

2.3730 


25.7000 

4.2392 

12.2000 

2.4064 


26.0000 

4 . 272 t 

12.4000 

2.4397 


26.3000 

4.3082 

12.6000 

2.4727 


26.6000 

4.3391 

12.8000 

2.5054 


26.9000 

4.3720 

13.0000 

2.5379 


27.2000 

4.4046 

13.2000 

2.5702 


27.5000 

4.4389 

13.4000 

2.6023 


27.8000 

4.4891 

13.6000 

2.6341 


28.1000 

4.5010 

13.8000 

2.6657 


28.4000 

4.5327 

14.0000 

2.6971 


28.7000 

4.5642 

14.2000 

2.7282 


29.0000 

4.5954 

14.4000 

2.7592 


29.3000 

4.6265 

14.6000 

2.7899 


29.6000 

4.6574 

14.0000 

2.8204 


29.9000 

4.6880 

15.0000 

2.8508 


30.2000 

4.7184 

15.2000 

2.8809 


30.5000 

4.7487 

15.4000 

2.9109 


30.8000 

4.7787 

15.6000 

2.9406 


31.1000 

4.8085 

15.8000 

2.9702 


31.4000 

4.8382 

16.0000 

2.9996 


31.7000 

4.8676 

16.2000 

3.0287 


32.0000 

4.8969 

16.4000 

3.0577 


32.3000 

4.9259 

16.6000 

3.0866 


32.6000 

4.9548 

16.8000 

3.1152 


32.9000 

4.9835 

17.0000 

3.1437 


33.2000 

5.0120 

17.2000 

3.1720 


33.5000 

5.0403 

17.4000 

3.2001 


33.8000 

5.0685 

17.6000 

3.2281 


34.1000 

5.0964 

17.8000 

3.2559 


34.4000 

5.1242 

18.0000 

3.2835 


34.7000 

5.1519 

18.2000 

3.3110 


35.0000 

5.1793 

18.4000 

3.3383 


35.3000 

5.2066 

18.6000 

3.3654 


35.6000 

5.2337 

18.8000 

3.3925 


35.9000 

5.2607 

19.0000 

3.4193 


36.2000 

5.2875 

19.2000 

3.4460 


36.5000 

5.3141 

19.4000 

3.4725 


36.8000 

5.3405 

19.6000 

3.4989 


37.1000 

5.3668 

19.8000 

3.5251 


37.4000 

5.3930 

20.0000 

3.5512 


37.7000 

5.4189 

20.3000 

3.5901 


38.0000 

5.4448 

20.6000 

3.6286 


38.3000 

5.4704 

20.9000 

3.6668 


38.6000 

5.4960 

21.2000 

3.7047 


38.9000 

5.5213 




Table 4. (continued) 


X , in . 

r , in . 


X , in . 

r , in . 

39.2000 

5.5465 


61.5000 

7,0755 

39.5000 

5.5716 


62.0000 

7.1032 

39.6000 

5.5965 


63.0000 

7.1581 

40.1000 

5.6213 


64.0000 

7.2119 

40.4000 

5.6459 


65.0000 

7.2649 

40.7000 

5.6704 


66.0000 

7.3169 

41.0000 

5.6947 


67.0000 

7.3681 

41.3000 

5.7189 


68.0000 

7.4183 

41.6000 

5.7430 


69.0000 

7.4677 

41.9000 

5.7669 


70.0000 

7.5163 

42.2000 

5.7907 


71.0000 

7.5640 

42.5000 

5.8144 


72.0000 

7.6108 

42.8000 

5.8379 


73.0000 

7.6569 

43.1000 

5.8613 


74.0000 

7.7021 

43.4000 

5.8845 


75.0000 

7.7465 

43.7000 

5.9076 


76.0000 

7.7902 

44.0000 

5.9306 


77.0000 

7.8330 

44.3000 

5.9535 


78.0000 

7.8752 

44.6000 

5.9762 


79.0000 

7.9165 

44.9000 

5.9988 


80.0000 

7.9571 

45.2000 

6.0212 


81.0000 

7.9970 

45.5000 

6.0436 


82.0000 

8.0362 

45.8000 

6.0658 


83.0000 

8.0747 

46.1000 

6.0879 


84.0000 

8.1124 

46.4000 

6.1099 


85.0000 

8.1495 

46.7000 

6.1317 


86.0000 

8.1858 

47.0000 

6.1535 


87.0000 

8.2215 

47.3000 

6.1751 


88.0000 

8.2565 

47.6000 

6.1966 


89.0000 

8.2908 

47.9000 

6.2180 


90.0000 

8.3245 

48.2000 

6.2392 


91.0000 

8.3575 

48.5000 

6.2604 


92.0000 

8.3899 

48.8000 

6.2814 


93.0000 

8.4216 

49.1000 

6.3023 


94.0000 

8.4526 

49.4000 

6.3231 


95.0000 

8.4830 

49.7000 

6.3438 


96.0000 

8.5128 

50.0000 

6.3644 


i 97.0000 

8. 5419 

50.5000 

6.3984 


' 98.0000 

8.5704 

51.0000 

6.4321 


99.0000 

8.5983 

51.5000 

6.4655 


100.0000 

8.6256 

52.0000 

6.4987 


101.0000 

8.6523 

52.5000 

6.5315 


102.0000 

8.6783 

53.0000 

6.5641 


103.0000 

8.7037 

53.5000 

6.5963 

j 

104.0000 

8.7286 

54.0000 

6.6283 


105.0000 

8.7528 

54.5000 

6.6600 


106.0000 

8.7764 

55.0000 

6.6914 


107.0000 

8.7993 

55.5000 

6.7225 


108.0000 

8.8217 

56.0000 

6.7533 


109.0000 

8.8435 

56.5000 

6.7839 


110.0000 

8.8646 

57.0000 

6.8143 


111.0000 

8.8850 

57.5000 

6.8443 


112.0000 

8.9049 

58.0000 

6.8741 


113.0000 

8.9240 

58.5000 

6.9036 


114.0000 

8.9424 

59.0000 

6.9329 


115.0000 

8.9601 

59.5000 

6.9619 


116.0000 

8.9768 

60.0000 

6.9907 


117.0000 

8.9924 

60.5000 

61.0000 

7.0192 

7.0475 


117.5410 

9.0000 




Table 5. Design Wall Coordinates for the Wj — 17 Nozzle 


X ) in . 

r , in . 

0.0000 

0.0607 

0.0018 

0.0607 

0.0066 

0.0607 

0.0135 

0.0607 

0.0235 

0.0607 

0.0379 

0.0608 

0.0585 

0.0608 

0.0883 

0.0609 

0.1311 

0.0610 

0.1928 

0.0611 

0.2817 

0.0613 

0.4096 

0.0615 

0.5934 

0.0620 

0.7934 

0.0629 

■ 0.9934 

0.0646 

1.1934 

0.0673 

1.3934 

0.0718 

1.5934 

0.0783 

1.7934 

0.0866 

1.9934 

0.0989 

2.1934 

0.1154 

2.3934 

0.1371 

2.5934 

0.1666 

2.7934 

0.2031 

2.9934 

0.2429 

3.1934 

0.2853 

3.3934 

0.3295 

3.5934 

0.3745 

3.7934 

0.4208 

3.9934 

0.4666 

4.1934 

0.5133 

4.3934 

0.5603 

4.5934 

0.6077 

4.7934 

0.6552 

4.9934 

0.7028 

5.1934 

0.7501 

5.3934 

0.7969 

5.5934 

0.8429 

5.7934 

0.8883 

5.9934 

0.9331 

6.1934 

0.9773 

6.3934 

1.0211 

6.5934 

1.0643 

6.7934 

1.1070 

6.9934 

1.1493 

7.1934 

1.1912 

7.3934 

1.2329 

7.5934 

1.2741 

7.7934 

1.3145 

7.9934 

1.3544 

8.2134 

1.3978 

8.5302 

1.4592 


X ) in . 

r , in . 

8.9302 

1.5355 

9.3302 

1.6102 

9.7302 

1.6836 

10.1302 

1.7555 

10.5302 

1.8266 

10.9302 

1.8960 

11.3302 

1 . 9645 ^ 

11.7302 

2.0322 

12.1302 

2.0903 

12.5302 

2 . 163 t - 

12.9302 

2 . 22 ti 

13.3302 

2.2920 

13.7302 

2.3847 

14.1302 

2.4160 

14.5302 

2.4791 

14.9302 

2.5407 

15.3302 

2.6006 

15.7302 

2.6599 

16.1302 

2.7189 

16.5302 

2.7773 

16.9302 

2.8346 

17.3302 

2.8912 

17.7302 

2.9473 

18.1302 

3.0027 

18.5302 

3.0576 

18.9302 

3.1118 

19.3302 

3.1655 

19.7302 

3.2188 

20.1302 

3.2715 

20.6582 

3.3403 

21.2582 

3.4175 

21.8582 

3.4936 

22.4582 

3.5686 

23.0582 

3.6426 

23.6582 

3.7155 

24.2582 

3.7874 

24.8582 

3.8585 

25.4582 

3.9286 

26.0582 

3.9980 

26.6582 

4.0665 

27.2582 

4.1341 

27.8582 

4.2010 

28.4582 

4.2671 

29.0582 

4.3325 

29.6582 

4.3971 

30.2582 

4.4611 

30.8582 

4.5243 

31.4582 

4.5868 

32.0582 

4.6488 

32.6582 

4.7099 

33.2582 

4.7706 

33.8582 

4.8307 



Table 5. (continued) 


75 


X, in. 

r, in. 


X, in . 

r , in . 

34.4582 

4.8900 


75.6182 

7.9191 

35.0582 

4.9488 


76.4162 

7.9622 

35.6582 

5.0069 


77.2182 

8.0050 

36.2582 

5.0646 


78.0182 

8.0474 

36.8582 

5.1217 


78.8182 

8.0894 

37.4582 

5.1781 


79.6182 

8.1311 

38.0582 

5.2341 


80.4182 

8.1724 

38.6582 

5.2895 


81.2182 

8.2134 

39.2582 

5.3444 


82.0182 

8.2540 

39.8582 

5.3987 


82.8182 

8.2942 

40.4582 

5.4526 


83.6182 

8.3341 

41.0582 

5.5059 


84.4182 

8.3736 

41.6582 

5.5586 


85.2182 

8 . 412 t 

42.2582 

5.6108 


86.0182 

8 . 4 S 1 S 

42.8582 

5.6626 


86.8182 

8 . 496 ft 

43.4582 

5.7141 


87.6182 

8 . 52 tft 

44.0582 

5.7649 


88.4182 

8 . 565 ft 

44.6582 

5.8152 


89.2182 

8.6032 

45.2582 

5.8650 


90.0182 

8.6402 

45.8582 

5.9144 


90.8182 

8 . 676 ft 

46.4582 

5.9633 


91.6182 

8.7131 

47.0582 

6.0120 


92.4182 

8.7490 

47.6582 

6.0599 


93.2182 

8 . 7846 . 

48.2582 

6.1076 


94.0182 

8.8198 

48.8582 

6.1556 


94.8182 

8.8546 

49.4582 

6.2026 


95.6182 

8.8891 

50.0582 

6.2495 


96.4182 

8.9233 

50.8182 

6.3080 


97.2182 

8.9570 

51.6182 

6.3689 


98.0182 

8.9904 

52.4182 

6.4292 


98.8182 

9.0235 

53.2182 

6.4890 


99.6182 

9.0562 

54.0182 

6.5480 


100.4180 

9.0885 

54.8182 

6.6069 


101.2180 

9.1205 

55.6182 

6.6641 


102.0180 

9.1521 

56.4182 

6.7210 


102.8180 

9.1833 

57.2182 

6.7775 


103.6180 

9.2142 

58.0182 

6.8331 


104.4180 

9.2447 

58.8182 

6.8884 


105.2180 

9.2749 

59.6182 

6.9430 


106.0180 

9.3047 

60.4182 

6.9969 


106.8180 

9.3342 

61.2182 

7.0504 


107.6180 

9.3633 

62.0182 

7.1032 


108.4180 

9.3920 

62.8182 

7.1556 


109.2180 

9.4204 

63.6182 

7.2072 


110.0180 

9.4484 

64.4182 

7.2584 


110.8180 

9.4760 

65.2182 

7.3090 


111.6180 

9.5033 

66.0182 

7.3591 


112.4180 

9.5303 

66.8182 

7.4086 


113.2180 

9.5568 

67.6182 

7.4576 


114.0180 

9.5831 

68.4182 

7.5061 


114.8180 

9.6089 

69.2182 

7.5541 


115.6180 

9.6344 

70.0182 

7.6016 


116.4180 

9.6595 

70.8182 

7.6486 


117.2180 

9.6843 

71.6182 

7.6952 


118.0180 

9.7087 

72.4182 

7.7413 


118.8180 

9.7328 

73.2182 

7.7867 


119.6180 

9.7565 

74.0182 

7.8318 


120.4180 

9.7799 

74.8182 

7.8756 


121.2180 

9.8028 


r 



Table 5. (continued) 


X, in. 

r, in. 

126.0180 

9.9332 

126.8180 

9.9536 

127.6180 

9.9737 

128.4180 

9.9935 

128.9060 

10.0053 


X, in. 

r, in. 

122.0180 

9.8254 

122.8180 

9.8477 

123.6180 

9.8696 

124.4180 

9.8912 

125.2180 

9.9123 


Table 0. Design Wall Coordinates for the CF^ Nossle 


X, in. 

r, in. 

0.0000 

0.2580 

0.1000 

0.2581 

0.2000 

0.2583 

0.3000 

0.2587 

0.4000 

0.2592 

0.5000 

0.2599 

0.6000 

0.2607 

0.7000 

0.2615 

0.8000 

0.2625 

0.9000 

0.2636 

1.0000 

0.2649 

1.1000 

0.2662 

1.3000 

0.2677 

1.3000 

0.2693 

1.4000 

0.2709 

1.5000 

0.2727 

1.6000 

0.2746 

1.7000 

0.2767 

1.8000 

0.2788 

1.9000 

0.2810 

2.0000 

0.2833 

2.1000 

0.2857 

2.2000 

0.2882 

2.3000 

0.2908 

2.4000 

0.2935 

2.5000 

0.2964 


X, in. 

r, in. 

2.6000 

0.2993 

2.7000 

0.3024 

2.8000 

0.3055 

2.9000 

0.3089 

3.0000 

0.3123 

3.1000 

0.3159 

3.2000 

0.3196 

3.3000 

0.3235 

3.4000 

0.3275 

3.5000 

0.3317 

3.6000 

0.3360 

3.7000 

0.3406 

3.8000 

0.3453 

3.9000 

0.3502 

4.0000 

0.3551 

4.1000 

0.3603 

4.2000 

0.3656 

4.3000 

0.3710 

4.4000 

0.3767 

4.5000 

0.3825 

4.6000 

0.3885 

4.7000 

0.3947 

4.8000 

0.4011 

4.9000 

0.4077 

5.0000 

0.4145 

5.1000 

0.4216 





Table 0. (continued) 


X , in . 

r , in . 


X , in . 

r , in . 

5.2000 

0.4288 


12.0000 

1.6670 

5.3000 

0.4362 


12.2000 

1.7173 

5.4000 

0.4439 


12.4000 

1.7679 

5.5000 

0.4519 


12.6000 

1.8186 

5.6000 

0.4600 


12.8000 

1.8694 

5.7000 

0.4684 


13.0000 

1.9205 

5.8000 

0.4771 


13.2000 

1.9716 

5.9000 

0.4861 


13.4000 

2.0229 

6.0000 

0.4953 


13.6000 

2.0743 

6.1000 

0.5048 


13.0000 

2.1258 

6.2000 

0.5146 


14.0000 

2.1774 

6.3000 

0.5247 


14.2000 

2.2290 

6.4000 

0.5352 


14.4000 

2.2807 

6.5000 

0.5459 


14.6000 

2.3324 

6.6000 

0.5570 


14.8000 

2.3842 

6.7000 

0.5685 


15.0000 

2.4360 

6.8000 

0.5803 


15.2000 

2.4877 

6.9000 

0.5925 


15.4000 

2.5394 

7.0000 

0.6051 


15.6000 

2.5912 

7.1000 

0.6181 


15.8000 

2.6431 

7.2000 

0.6316 


16.0000 

2.6950 

7.3000 

0.6455 


16.2000 

2.7469 

7.4000 

0.6599 


16.4000 

2.7980 

7.5000 

0.6749 


16.6000 

2.0506 

7.6000 

0.6903 


16.8000 

2.9022 

7.7000 

0.7063 


17.0000 

2.9535 

7.8000 

0.7228 


17.2000 

3.0046 

7.9000 

0.7398 


17.4000 

3.0554 

8.0000 

0.7573 


17.6000 

3.1060 

8.1000 

0.7752 


17.8000 

3.1563 

8.2000 

0.7936 


18.0000 

3.2063 

8.3000 

0.8124 


18.2000 

3.2561 

8.4000 

0.8316 


18.4000 

3.3056 

8.5000 

0.8511 


18.6000 

3.3548 

8.6000 

0.8710 


18.8000 

3.4037 

8.7000 

0.8913 


19.0000 

3.4522 

8.8000 

0.9119 


19.2000 

3.5006 

8.9000 

0.9327 


19.4000 

3.5486 

9.0000 

0.9539 


19.6000 

3.5963 

9.1000 

0.9753 


19.8000 

3.6437 

9.2000 

0.9970 


20.0000 

3.6909 

9.3000 

1.0189 


20.3000 

3.7610 

9.4000 

1.0411 


20.6000 

3.8305 

9.5000 

1.0635 


20.9000 

3.8993 

9.6000 

1.0861 


21.2000 

3.9675 

9.7000 

1.1088 


21.5000 

4.0350 

9.8000 

1.1318 


21.8000 

4.1019 

9.9000 

1.1549 


22.1000 

4.1682 

10.0000 

1.1782 


22.4000 

4.2338 

10.2000 

1.2252 


22.7000 

4.2988 

10.4000 

1.2727 


23.0000 

4.3631 

10.6000 

1.3208 


23.3000 

4.4269 

10.8000 

1.3693 


23.6000 

4.4900 

11.0000 

1.4181 


23.9000 

4.5526 

11.2000 

1.4673 


24.2000 

4.6144 

11.4000 

1.5168 


24.5000 

4.6757 

11.6000 

1.5668 


24.8000 

4.7365 

11.8000 

1.6168 


25.1000 

4.7966 




Table 6. (continued) 


X, in. 

r, in. 


X, in. 

r, in. 

25.4000 

4.8561 


48.0004 

8.0237 

25.7000 

4.9152 


48.5004 

8.0705 

26.0000 

4.9735 


49.0004 

8.1166 

26.3000 

5.0314 


49.5004 

8.1619 

26.6000 

5.0886 


50.0004 

8.2064 

26.9000 

5.1453 


51.0004 

8.2931 

27.2000 

5.2015 


52.0004 

8.3769 

27.4999 

5.2572 


53.0004 

8.4577 

27.7999 

5.3123 


54.0004 

8.5358 

28.0999 

5.3668 


55.0004 

8.6110 

28.3999 

5.4208 


56.0008 

8.6836 

28.6999 

5.4743 


57.0008 

8.7535 

28.8999 

5.5097 


58.0008 

8.8210 

29.1999 

5.5623 


59.0008 

8.8860 

29.4999 

5.6145 


60.0008 

8.9485 

29.7999 

5.6661 


61.0008 

9.0087 

29.9999 

5.7002 


62.0008 

9.0665 

30.4999 

5.7846 


63.0008 

9.1222 

30.9999 

5.8677 


64.0012 

9.1756 

31.4999 

5.9494 


65.0012 

9.2267 

31.9999 

6.0298 


66.0012 

9.2758 

32.4999 

6.1091 


67.0012 

9.3227 

32.9999 

6.1870 


68.0012 

9.3678 

33.4999 

6.2637 


69.0012 

9.4108 

33.9999 

6.3393 


70.0012 

9.4519 

34.4999 

6.4136 


71.0012 

9.4912 

34.9999 

6.4868 


72.0012 

9.5287 

35.4999 

6.5588 


73.0016 

9.5643 

35.9999 

6.6297 


74.0016 

9.5983 

36.4999 

6.6994 


75.0016 

9.6306 

36.9999 

6.7679 


76.0016 

9.6612 

37.4999 

6.8354 


77.0016 

9.6903 

37.9999 

6.9017 


78.0016 

9.7177 

38.4999 

6.9670 


79.0016 

9.7437 

38.9999 

7.0313 


80.0016 

9.7682 

39.5000 

7.0945 


81.0016 

9.7913 

40.0000 

7.1567 


82.0016 

9.8130 

40.5000 

7.2178 


83.0020 

9.8334 

41.0000 

7.2780 


84.0020 

9.8524 

41.5000 

7.3372 


85.0020 

9.8701 

42.0000 

7.3954 


86.0020 

9.8867 

42.5000 

7.4527 


87.0020 

9.9019 

43.0000 

7.5091 


88.0020 

9.9157 

43.5000 

7.5645 


89.0020 

9.9283 

44.0000 

7.6190 


90.0020 

9.9396 

44.5000 

7.6726 


91.0024 

9.9501 

45.0000 

7.7253 


92.0024 

9.9597 

45.5000 

7.7772 


93.0024 

9.9685 

46.0000 

7.8282 


94.0024 

9.9767 

46.5004 

7.8783 


95.0024 

9.9845 

47.0004 

7.9276 


95.5024 

9.9882 

47.5004 

7.9761 


95.7024 

9.9897 
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Fig. 1.1 Pitot'pressure profiles in test section of Nitrogen IHumd. 
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Fig. i.a Pitot-pressuxe profiles in test section of CF^ Tunnel. 
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a) physical domain 



b) computational domain 


Fig. 3.0 Grid transformation. 
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Pig* 0.0 Design wall temperature distributions. 
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Fig. 0.1.1 Air nozzle contour compared to existing Mach 10 nozzle. 
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Fig. 6*1.8 Total pressure profile downstream of normal shock in exit plane 

I of Air nosale. 

1 

I 






















Fig. 6.2.0 Navier-Stokes Mach number odn<min for JV, - 17 nozzle. 





Fig. 6.3.12 Centerline pressure ratio for jVj — 17 nozzle (second wall temp.). 
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Fig. 6.3 .13 Movement of inflection point due to discretization. 







